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pH . ABSTRACT 

I We report a new analysis of the stellar dynamics in the Galactic centre, based on 

. improved sky and line-of-sight velocities for more than one hundred stars in the central 

' few arcseconds from the black hole candidate SgrA*. The main results are: 

O^l i • Overall the stellar motions do not deviate strongly from isotropy. For those 

32 stars with a determination of all three velocity components the absolute, line of 
I sight and sky velocities are in good agreement, consistent with a spherical star cluster. 

, Likewise the sky-projected radial and tangential velocities of all 104 proper motion 

' stars in our sample are also consistent with overall isotropy. 

I • However, the sky-projected velocity components of the young, early type stars 

I . in our sample indicate significant deviations from isotropy, with a strong radial de- 

' pendence. Most of the bright Hel emission line stars at separations from 1" to 10" 

', from SgrA* are on tangential orbits. This tangential anisotropy of the Hel stars and 

■ most of the brighter members of the IRS 16 complex is largely caused by a clockwise 

(on the sky) and counter-rotating (line of sight, compared to the Galaxy), coherent 
rotation pattern. The overall rotation of the young star cluster probably is a remnant 
of the original angular momentum pattern in the interstellar cloud from which these 

O I stars were formed. 

^ . • The fainter, fast moving stars within « 1" from SgrA* appear to be largely 

^ ' moving on radial or very elliptical orbits. We have so far not detected deviations from 

. . I linear motion (i.e. acceleration) for any of them. Most of the SgrA* cluster members 

■ also are on clockwise orbits. Spectroscopy indicates that they are early type stars. 
' We propose that the SgrA* cluster stars are those members of the early type cluster 
, that happen to have small angular momentum and thus can plunge to the immediate 

■ vicinity of SgrA*. 

• We derive an anisotropy-independent estimate of the Sun-Galactic centre dis- 
tance between 7.8 and 8.2 kpc, with a formal statistical uncertainty of ±0.9 kpc. 

• We explicitly include velocity anisotropy in estimating the central mass distribu- 
tion. We show how Leonard-Merritt and Bahcall-Tremaine mass estimates give system- 
atic offsets in the inferred mass of the central object when applied to finite concentric 
rings for power law clusters. Corrected Leonard-Merritt projected mass estimators and 
Jeans equation modelling confirm previous conclusions (from isotropic models) that a 
compact central mass concentration (central density >10 ^^'^Mq pc is present and 
dominates the potential between 0.01 and 1 pc. Depending on the modelling method 
used the derived central mass ranges between 2.6 and 3.3 x 10 ^Mq for Rq = 8.0 kpc. 



1 INTRODUCTION 

High spatial resolution observations of the motions of gas and stars have in the past decade substantially strengthened the 
evidence that central dark mass concentrations reside in many and perhaps most nuclei of nearby galaxies (Kormendy and 
Richstone 1995, Magorrian et al. 1998, Richstone et al. 1998). These dark central masses are very likely massive black holes. 
The most compelling evidence for this assertion comes from the dynamics of water vapor maser cloudlets in the nucleus of 
NGC 4258, and from the stellar dynamics in the centre of our own Galaxy (Greenhill et al. 1995, Myoshi et al. 1995, Eckart 
and Genzel 1996, 1997, Genzel et al. 1997, Ghez et al. 1998). In both cases the gas and stellar dynamics indicate the presence 
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Figure 1. Examples of the proper motions derived from the 1992 to 1998 data NTT data sets obtained with the MPE SHARP camera. 
Shown are position-time plots for source SI close to SgrA* (left) and the bright Hel emission line star IRS 16 C (right). The two panels 
show x-(=RA) and y-(=Dec) position offsets, along with the best-fitting proper motions (straight lines). For e3x:h epoch the average and 
1(7 dispersion of a number position measurements are plotted. 



of an unresolved central mass whose density is so large that it cannot be stable for any reasonable length of time unless it is 
in form of a massive black hole (Maoz 1998). 

The case of the Galactic centre is particularly intriguing, as it is very close (8 kpc). With the highest spatial resolution 
observations presently available in the near-infrared (<0.1"), spatial scales of a few light days can bo probed. Measurements 
of both line-of-sight velocities (through Dopplcr shifts in spectral lines) and sky/proper motions arc available and pose very 
strong constraints on the central mass concentration. The following results have emerged. 

• The mean stellar velocities (or velocity dispersions) follow a Kepler law ((v ^) oc R ~^) from projected radii R from 
0.1" to Ri20", providing compelling qualitative evidence for the presence of a central point mass (Sellgren et al. 1990, Krabbe 
et al. 1995, Hallor ct al. 1996, Genzel ct al. 1996, 1997, Eckart and Gonzcl 1996, 1997, Ghcz ct al. 1998). 

• The positions of the dynamical centre (maximum velocity dispersion) and of the maximum stellar surface number 
density agree with the position of the compact radio source SgrA* (size less than a few AU, Lo et al. 1998, Bower and Backer 
1998) to within ±0.1" (Ghez et al. 1998). 

• Projected mass estimators and Jeans equation modeling of the stellar velocity data indicate that the central mass ranges 
between 2.2 and 3 x 10 ''Mq. It has a mass to luminosity ratio of M/L > WQMq / Lq and a density > 2 x 10 ^^Mq pc ~^ 
(Genzel et al. 1997, Ghez et al. 1998). 

The Galactic centre mass modeling has so far assumed that the stellar velocity ellipsoid is isotropic. An initial compari- 
son of line-of-sight and proper motion velocity dispersions indeed suggests that there are no coarse deviations from isotropy 
(Eckart and Genzel 1996, 1997). However, to make a more detailed assessment, it is necessary to obtain more accurate stellar 
motions than were available two years ago. These improved motions -for line-of-sight and sky components - are now in hand 
and will be analysed in the present paper. 



2 OBSERVATIONS AND DATA ANALYSIS 
2.1 Proper Motions 

In our earlier papers (Eckart et al. 1992, 1993, 1994, 1995) we have described the data acquisition and reduction that allowed 
us to obtain stellar positions with a precision of «10 milli-arcseconds per measurement (Eckart and Genzel 1996, 1997, Genzel 
et al. 1997). We used the MPE-SHARP camera (Hofmann et al. 1993) on the 3.5m New Technology Telescope (NTT) of the 
European Southern Observatory (ESO). SHARP contains a 256^ pixel NICMOS 3 detector. Each pixel projects to 25 or 50 
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Figure 2. K-band spectra of early type stars in the central f» 10" obtained in spring 1996 with the MPE 3D spectrometer on the 
ESO-MPG 2.2m telescope . The spectral resolving power is A/AA = 2000, Nyquist sampled at twice that resolution. Most of the spectra 
were obtained by subtracting from the 'on-star' spectrum (0.3" to 1.2" aperture) an 'off-star' spectrum scaled to the same area, in 
order to remove local nebular emission. The name is given for each star, as is the stellaj: velocity in km/s, with the 1 a uncertainty in 
parentheses. 



milli-arcseconds on the sky in order to (over-)sample the « 0.15" FWHM diffraction limited image of the NTT in K-band. 
The raw data for each data set consist of several thousands of short exposure frames (0.3 to 0.5 seconds integration time). 

First wo processed the data from nights with very good seeing (0.4 to 0.8" ) in the standard manner (dead pixel correction, 
sky subtraction, flat-fielding etc.). Next we co-added with the simplo-shift-and-add algorithm (SSA: for details see Christou 
1991 and Eckart et al. 1994). The individual short exposure frames typically contain only a few bright speckles so that the 
SSA algorithm is well suited for our purpose. The bright infrared sources IRS7 or IRS16NE serve as reference sources. For the 
present study we analysed «82 independent data sets from a total of 9 observing runs in 1992.25, 1992.65, 1993.65, 1994.27, 
1995.6, 1996.25, 1996.43, 1997.54 and 1998.37. For the central arcsccond region around SgrA* (the 'SgrA* cluster') wo also 
analysed an additional data set from 1999.42. Eckart and Genzel (1996, 1997) and Genzel et al. (1997 have previously analysed 
and discussed the data until 1996.43. 



The diffraction limited core of the stellar SSA images contains up to 20 percent of the light. Determinations of the 
relative pixel offsets from IRS16NE in raw SSA images or from diffraction limited maps after removing the seeing halo give 
consistent results (see Eckart and Genzel 1996, 1997 for details). Those 'CLEANod' SSA maps produce similar results as 
other data reductions (Knox-Thompson, triple correlation) but give much higher dynamical range (see Eckart et al. 1994 for 
a detailed discusssion). This is essential in the crowded Galactic centre region where the difference between the brightest and 
faintest stars in our images is almost 10 magnitudes. We solved for the relative offsets, rotation angle, and for the linear and 
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Figure 3. see Fig. (g) 

quadratic distortions between individual frames from an over-determined system of non-linear equations for a reference list of 
relatively isolated bright stars (Eckart and Genzel 1996, 1997). Our final near- infrared reference frame is tied to an accuracy 
of ±30 milli-arcsecs to the VLA radio frame through 5 stars that show SiO and H 2O maser activity and are common to both 
wavelength bands (Menten et al. 1997). The resulting combined systematic errors in our proper motion velocity estimates 
probably are about 30 km/s. In the best cases these systematic effects now dominate the error budget. 

The new 1997 and 1998 data sets are in excellent agreement with the extrapolation of the data we have published before 
and significantly improve the uncertainties. As examples we show in Fig. (|l|) the relative RA- and Dec-position offsets as a 
function of time between 1992 and 1998 for two selected stars. IRS16 C is a bright and isolated, Hel emission line star (Krabbe 
et al. 1991, 1995). Its position vs. time diagram in Fig. (0) is an example of the quality of the data on bright isolated stars. 
SI is a faint star in the 'SgrA* cluster' that is very close to SgrA* (~0.1"). It shows the fastest proper motion (~1470 km/s) 
in the entire sample. 



2.2 3D Spectroscopy 

We observed the Galactic centre with the MPE-3D near-infrared spectrometer (Weitzel et al. 1996) in conjunction with 
the tip-tilt adaptive optics module ROGUE (Thatte et al. 1995). 3D is a field imaging spectrometer which obtains spectra 
simultaneously for 256 spatial pixels covering a square region of the sky (16x16 pixels). The fill factor is over 95For further 
details of the instrument and data reduction we refer to Weitzel et al. (1996). We observed the Galactic centre in March 1996 
at the 2.2m ESO-MPG telescope on La Silla, Chile. During the run the seeing on the seeing monitor ranged between 0.3" and 
0.8". The pixel scale was 0.3". We observed the short-wavelength part of the K-band (1.9-2.2/im) at A/AA=2000, Nyquist 



© 0000 RAS, MNRAS 000, 000-000 



Stellar Dynamics in the Galactic centreiProper Motions and Anisotropy 5 



sampled with two settings of a piezo driven flat mirror. We covered the central ~16"xl0" centered on SgrA* by an overlapping 
set of frames, each with a field of view of 4.8"x4.8". At each position we set up a sequence on-source (piezo step 1), off-source 
(piezo step 1), on-source (piezo step 2), off-source (piezo step 2) etc. with an integration time per step of 200 seconds. Due to 
the combined effects of seeing and pixel scale the resulting FWIfM spatial resolution of the final combined data set was 0.6". 
We employed the standard 3D data analysis package (based on GIPSY: van der Hulst et al. 1992). We performed wavelength 
calibration, sky subtraction, spectral and spatial flat flelding, dead and hot pixel correction and division by a reference stellar 
spectrum obtained during the observations. We corrected for the effects of a spatially varying fringing or 'channel' spectrum 
due to interference in the saphire coating of the NICMOS 3 detector by applying suitable flat fields from a set of flatflelds 
at different settings of the piezo mirror. Based on obervations of calibration lamps and OH sky lines during the different 
observing nights the flnal velocity calibration is accurate to ±16fcm/s. 

In the observed field we identified 21 emission line stars from continuum subtracted images of the 2.058/im n=2 n=2 
and 2.11/jm n=4 ^'^S - n=3 ^'^P Hel lines and the 2.166^m n=7-4 HI (Br7) line. Most of the stars are identical with those 
found in the «1", R=1000 3D data set of Krabbe et al. (1995) and Genzel et al. (1996) ( see also Blum, Sellgren and dePoy 
1995 a,b and Haller et al.l986) but the resolution, quality and nebular rejection is now much superior. Three new stars were 
identified: (-2.1", -4.1") (13S SE), {+1.6" , +0.3") (16 CC) and (-8.3", -5.7") (aU offsets are in RA and Dec from SgrA*). 
We extracted from the data cube spectra of individual stars by typically averaging 3 to 16 pixels on the star, for effective 
apertures between 0.3" and 1.2". In most cases, we subtracted a suitable 'off-star' spectrum (scaled to the same aperture 
area) to eliminate the effect of local nebular line emission. Fig. (^ and Fig. (^ show the final spectra for 18 of the 21 stars. 

To determine stellar velocities we fitted Gaussians to the 2.058^m and 2.11^m Hel lines and the 2.166/im Br7 line. In a few 
cases we also used in addition the 2.137/2.143^m Mgil lines and the 2.189/im Hell line. A number of the stars clearly display 
P-Cyg profiles in the 2.058/im Hel transition (Fig. (g) and Fig. ^) . In these cases we fitted the profile with a double Gaussian 
(absorption and emission). As the absorption structure is well resolved an unambiguous emission line centroid (assumed to be 
the stellar velocity) can thus be easily obtained. For most stars we determined the final stellar velocities from averaging the 
values obtained from 3 (or 4) lines. The agreement between the fits to the different lines is generally good or even excellent. 
We list the final velocities in Table 1 and the insets of Fig. (^ and Fig. (^). The new velocity determinations agree with those 
of Genzel et al. (1996) but the uncertainties are typically half of those in our earlier work. The best cases have an uncertainty 
of ±25 km/s. 



2.3 A Homogenized Data Set 

To obtain a homogenized 'best' data set of stellar velocities for further analysis we combined the new 1992-1998 NTT proper 
motions and 2.2m line-of-sight velocity data described in the last two paragraphs with the 1995-1997 Keck proper motion data 
of Ghez et al. (1998) and with other relevant line-of-sight velocity data sets (see Genzel et al. 1996 and references therein). 
Table 1 gives the results. The following explanations and comments for Table 1 are in order. 

• Columns 1 through 3 contain the projected separation R, x(=RA)-offset and y (=Dec)-offset between star and SgrA* 
(epoch 1994/1995, all in arcseconds). As for the NTT, the Keck astrometry is established with H20/SiO maser stars that are 
visible in both wavelength bands and is accurate to ±10 milli-arcseconds (see Ghez et al. 1998, Menten et al. 1997). 

• Columns 4 through 15 contain the x- and y-proper motions and their respective la errors (km/s, for a Sun-Galactic centre 
distance of Rq = S.Okpc). Whenever two measurements are available, we first list the velocity from the Keck observations 
and then that from the NTT observations. Columns 12 through 15 give the final combined proper motions obtained from 
averaging motions from the two sets (if available) with 1/(t^ weighting. The agreement between Keck and NTT data sets 
generally is very good (see discussion by Ghez et al. 1998) and is in accordance with the measurement uncertainties. We thus 
assume that the final measurement error of the combined set is given by 1/ ^y (l/c"i ± l/o"!)- We have found, however, from 
a comparison of the two data sets that stars with >200 km/s velocity uncertainty at _R > 1" (and >400 km/s at 7? < 1") are 
highly unreliable. We therefore have eliminated such stars from the final set. 

• Columns 16 and 17 give the line of sight velocity and its la error. 

• Column 18 assigns a weight to each data point, based on its reliability, the velocity errors and the agreement between 
different data sets. The weight is approximately proportional to 1/error^ as appropriate for white noise but 'quantized' 
in order to not place too much weight on the few data points with the smallest statistical errors. While this weighting 
scheme is subjective it is in our opinion a fair representation of the quality of the different data points in the presence of 
significant systematic errors. We have also tried another, more formal weighting scheme. Here we have assigned the weight 
w = 1/(1 ± (error/o"o)^), where 'error' is the x/y-averaged proper motion velocity uncertainty, ao = 250f?~°'^ is a measure 
of the sample dispersion at R. This weighting scheme also gives essentially the same weights for different data points as long 
as their individual errors are much smaller than the sample velocity dispersion. For large errors the weight scales as 1/error^ 
as for white noise. We have applied both weighting schemes in the various estimates discussed in the text. The results are 
basically identical. 

• Column 19 lists the 'popular' name of the star. 

• Column 20 lists an identification; 'p' stands for a star with a measured proper motion, 'early' denotes that the star is 
a young early type star (e.g. Hel/HI emission lines), and 'CO' denotes that the star is a late type star. 
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265 


33 


94 


46 


-60 


30 




16NW' 


p/early 


all 


9.86 





1.20 


-1.20 


-0.11 


120 


170 


330 


60 


-125 


200 


365 


61 


17 


130 


347 


43 








W8 


P 


K/N 


12.52 





1.29 


1.22 


0.44 


-370 


60 


380 


40 


-301 


52 


280 


67 


-330 


39 


353 


34 


160 


25 


1 


16C 


p/early 


aU 


9.55 





132 


-0.85 


-1.00 


220 


100 


20 


100 


259 


150 


-260 


120 


232 


83 


-95 


77 


-1 


30 


1 


W13 


pCO 


aU 


13.3 




1J4 


-0.84 


-1.05 


232 


115 


-287 


77 










232 


115 


-287 


77 






0.5 


w217 


P 


NTT 






1J6 


-1.03 


-0.88 


-410 


120 


-260 


140 


-494 


130 


-410 


68 


-449 


88 


-382 


61 






1 


W12 


P 


K/N 


13.8 




137 


-1.34 


-0^9 


-20 


70 


-100 


60 


-85 


120 


-272 


84 


-37 


60 


-158 


49 


270 


70 


1 


WIO 


p/early 


aU 


12.48 





141 


-1.31 


0.52 


30 


160 


-30 


90 


-105 


118 


-193 


86 


-58 


95 


-116 


62 






1 


W4 


P 


K/N 


13.9 




1.44 


1.06 


-0.98 


170 


60 


150 


40 


360 


70 


70 


70 


250 


46 


130 


35 


450 


30 


1 


16SW 


p/early 


aU 


9.61 


2 (doul 
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o 
o 
o 

> 

CO 



> 

CO 

o 
o 
p 

o 
o 

? 

o 
o 
o 



Table 1 (cont'd 



(SgrA*) 


Ax(SgrA*) 


Ay(SgrA*) 


Vx 






Avy 




9 


10 


11 

Av, 


Vrf 


AVrf 


Vyf 


AVy, 






weight 


name 


type 


source 


K-mag 


variabililty 




(arsccs) 














all velocities are in km/s 


























1^5 


0.55 


-145 


-200 


50 


-80 


50 










-200 


50 


-80 


50 






1 






Keck 


12 




1.63 


038 


-1.58 


20 


50 


-60 


120 










20 


50 


-60 


120 






0.5 






Keck 


134 




1.64 


0.41 


1.59 


240 


80 


160 


50 


299 


77 


29 


53 


271 


55 


98 


36 






1 






K/N 


12.27 





1.65 


-1.65 


0.14 


160 


120 


-200 


60 


-94 


126 


-185 


95 


39 


120 


-196 


51 






1 


W7 




K/N 


13.3 




1.67 


-1.65 


0.28 


-107 


114 


-127 


86 










-107 


114 


-127 


86 






0.5 






NTT 


12.86 





1.70 


-1.63 


-0.50 


70 


60 


-5 


80 


150 


135 


-83 


77 


83 


55 


-45 


55 








W14 




K/N 


12.52 





1.72 


-0.90 


-1.46 


-20 


80 


10 


90 


14 


69 


-30 


50 





52 


-21 


44 






J 






Keck 


11.42 





1.73 


0.76 


-1.55 


-160 


40 


-210 


70 


71 


70 


-394 


41 


-103 


35 


-348 


35 


-209 


30 


1 




pCO 


K/N 


11.14 





1.81 


1.69 


-0.66 


-70 


90 


240 


50 


132 


241 


-129 


85 


-45 


84 


145 


43 


-175 


30 


1 




pCO 


K/N 


13.4 




1.89 


0.63 


-1.78 


1 10 


50 


-260 


50 










110 


50 


-260 


50 






1 






NTT 


11.5 





2.02 


-1.80 


0.93 


-40 


50 


50 


50 


12 


112 


-2 


85 


-31 


46 


37 


43 


-93 


20 


1 


29S 


pCO 


all 


10.95 





2.07 


-0.54 


2.00 


140 


140 


200 


40 










140 


140 


200 


40 






0.5 






Keck 


13.8 




2.09 


1.46 


-1.49 


140 


70 


230 


40 










140 


70 


230 


40 






1 






Keck 


11.67 





2.09 


1.91 


-0.86 


240 


110 


80 


80 


36 


128 


-17 


53 


154 


84 


13 


44 






1 






K/N 


12.9 




2.10 


1.60 


-1.36 


200 


60 


230 


50 










200 


60 


230 


50 












Keck 


11.67 





2.10 


1.06 


1.81 


-60 


90 


220 


80 


-88 


64 


42 


80 


-79 


52 


132 


57 






J 




p 


K/N 


13.2 




2.11 


-1.94 


0.82 


150 


110 


120 


90 





160 


243 


117 


102 


91 


166 


71 








W2 




K/N 


11.9 




2.12 


2.06 


0.51 


-120 


50 


170 


60 


5 


50 


250 


30 


-58 


35 


234 


27 


210 


70 


1 


16CC 


p/earlv 


all 


10.15 





2.13 


-1.59 


1.41 


200 


100 


-140 


100 


174 


94 


-58 


75 


186 


68 


-88 


60 


-110 


150 


1 


29N 


p/earlv 


all 


9.74 





2.13 


2.03 


-0.63 


-130 


80 


110 


70 










-130 


80 


110 


70 






1 






Keck 


11.5 




2.16 


-0.02 


-2.16 


-90 


140 


-100 


90 


92 


51 


-259 


30 


70 


48 


-243 


28 








33 N 


p/early? 


Keck 


10.84 





2.20 


1.86 


-1.16 


170 


60 


240 


40 


278 


31 


50 


28 


255 


28 


113 


23 


420 


60 


J 


16SE1 


p/earlv 


all 


10.56 





2.23 


-0.91 


2.04 


-200 


90 


10 


40 


-310 


50 





50 


-284 


44 


6 


31 


-180 


100 




29NE1 


p/earlv 


all 


11.59 


1 


2.24 


-1.78 


-136 


-342 


133 


-63 


92 










-342 


133 


-63 


92 






0.25 


W15 




NTT 






2.26 


1.27 


-1.87 


260 


60 


180 


40 


383 


48 


-138 


40 


335 


38 


21 


28 






1 




p/early? 


K/N 


10.44 





232 


-0.90 


-2.14 


-320 


230 


140 


140 


-208 


62 


-68 


40 


-216 


60 


-52 


38 


102 


25 


1 




pCO 


all 


12.69 





233 


0.53 


2.27 


-180 


180 


90 


40 


-258 


102 


39 


43 


-239 


88 


66 


29 


-107 


40 


1 




pCO 


all 


1232 





236 


-1.70 


-1.65 


420 


110 


80 


100 


-8 


124 


-58 


60 


231 


82 


-22 


51 






1 






K/N 


13.2 




239 


2.37 


-0^9 


-140 


60 


330 


40 


-42 


106 


346 


234 


-116 


52 


330 


39 


-143 


30 


1 




pCO 


aU 


12.53 





246 


-2.28 


-0.93 


-140 


220 


110 


230 


-186 


170 


166 


208 


-169 


134 


141 


154 






0.25 


W16 


P 


K/N 


13.5 




2.59 


0.79 


-246 


-200 


290 


200 


200 


437 


229 


-8 


100 


193 


180 


33 


90 






0.25 




P 


K/N 


13.7 




2.60 


0.73 


2.50 


340 


120 


-420 


100 


442 


499 


-825 


489 


346 


117 


-436 


98 


-100 


30 


0.25 




pCO 


aU 


1238 





2.76 


0.27 


2.74 


-170 


270 


-690 


190 


-385 


82 


-21 


78 


-366 


79 


-117 


72 






0.5 




P 


K/N 


12.8 





2.77 


-1.80 


-2.11 


-192 


49 


34 


40 










-192 


49 


34 


40 






1 




P 


NTT 


12.89 





2.88 


-1.14 


2.65 


-130 


50 


-140 


30 










-130 


50 


-140 


30 


17 


30 


1 




pCO 


NTT/LS 


11.44 


1 


2.92 


2.91 


-OOO 


-360 


60 


140 


50 


-197 


102 


28 


26 


-318 


52 


52 


23 






1 




p 


K/N 


12.55 





2.98 


1.15 


2.75 


120 


210 


260 


110 


165 


65 


55 


40 


161 


62 


79 


38 






1 




p 


K/N 


11.89 





3.07 


-1.13 


2.85 


-184 


33 


-196 


57 










-184 


33 


-196 


57 






0.5 




p 


NTT 


s 




3.09 


2.89 


1.10 


160 


90 


-290 


30 


242 


95 


-267 


31 


199 


65 


-279 


21 


-15 


25 


1 


16NE 


p/early 


aU 


8.76 





3.12 


1.28 


2.84 


103 


48 


118 


30 










103 


48 


118 


30 






0.5 




P 


NTT 


X 





(1) 

O 
3 

o 

cc 



cr 
Q 
o 



O 

(1) 

C3 

CD 

t 

-i 
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o 
S 

Co 

a 

Si- 

Co 

O 

I 



Table 1 (cont'd 



@ 

O 
O 

o 
o 

> 

CO 



> 

CO 

o 
o 
o 

o 
o 
o 

i 

o 
o 



1 

[SgrA*) 


2 

Ax(SgrA*) 

(arsccs) 


3 

Ay(SgrA*) 


4 

Vx 


5 




7 

Avy 


g 

Vx 


9 10 11 

AVx Vy AVy 

all Tclocitics are in km/s 


12 
Vrf 


13 

AVrf 


3.12 


3.07 


0.56 


210 


90 


80 


70 


394 


209 118 142 


239 


83 


3.21 


2.97 


-1.20 


50 


90 


260 


130 


192 


62 79 61 


147 


51 


3.26 


3.26 


0.08 


90 


100 


-280 


70 






90 


100 


3.28 


2.01 


-2.59 


-90 


260 


200 


240 


-30 


148 -193 124 


-44 


129 


3.34 


0.50 


-3.30 


211 


60 


-60 


30 






211 


60 


3.37 


-2.85 


-1.80 


















3.45 


1.63 


-3.04 


1 10 


158 


-21 


48 






110 


158 


3.46 


-1.80 


-2.95 


213 


60 


172 


35 






213 


60 


3.47 


-0.50 


-3.43 


60 


60 


30 


35 






60 


60 


3.51 


3.30 


-1.20 


-80 


70 


20 


60 






-80 


70 


3.59 


3.40 


-1.16 


-60 


70 


-10 


120 






-60 


70 


3.59 


3.03 


-1.92 


-120 


200 


240 


130 






-120 


200 


3.65 


0.79 


3.56 


113 


62 


-94 


51 






113 


62 


3.70 


2.48 


-2.74 


-270 


170 


170 


110 


-140 


70 50 40 


-159 


65 


3.80 


-3.62 


-1.15 


8 


182 


583 


149 






8 


200 


4.08 


2.55 


3.19 


309 


159 


429 


105 






309 


159 


4.14 


3.78 


1.70 


-160 


100 


-360 


80 






-160 


100 


4.32 


3.65 


-2.32 


-100 


130 


-1 10 


170 






-100 


130 


4.39 


-4.08 


1.62 


-80 


50 


-130 


40 






-80 


50 


4.44 


-2.31 


3.79 


170 


40 


115 


45 






170 


40 


4.56 


-2.10 


-4.05 


















4.67 


4.65 


0.45 


















4.69 


-3.70 


2.88 


-60 


100 


-200 


80 






-60 


100 


4.81 


3.50 


-330 


















4.89 


-4.20 


-2.50 


















4.98 


-4.60 


-1.90 


















5.04 


5.00 


-0.60 


















5.16 


2,33 


4.60 


371 


156 


225 


69 






371 


156 


5J0 


5.30 


0.10 


















5.44 


1.90 


-5.10 


















5.57 


5.55 


0.45 


















558 


-1.40 


-540 


















5.75 


0.22 


5.75 


100 


67 


-118 


35 






100 


67 


5.77 


-0.90 


-5.70 


















5.77 


-5.70 


0.90 


















5.79 


3.00 


-4.95 


















5.88 


3.10 


5.00 


















5.95 


3.90 


4.50 


















6.04 


-5.30 


2.90 


















6.04 


-3.80 


-4.70 



















14 


15 


16 


17 


18 


19 


20 


21 


22 


23 








Vj 


Av. 


weight 


name 


type 


source 


K-mag 


variabililty 




87 


63 






1 




p/early? 


K/N 


11.93 







111 


55 


230 


90 


1 


16SE2 


p/early 


all 


11.75 





hch 


-280 


70 






0.5 




p/early? 


Keck 


12.43 





-111 


110 






0.25 




p 


K/N 


13.7 




on, 


-60 


30 


160 


60 




33E 


p/early 


NTT/LS 


9.9 









10 


60 


1 


13E 


early 


LaSilla96 


10.26 






-21 


48 






0.25 




P 


NTT 


X 






172 


35 






1 




p 


NTT 


11.71 





30 


35 


82 


25 


1 


33W 


pCO 


NTT/LS 


10.81 





'kart. 


20 


60 






0.5 




P 


Keck 


11.48 





-10 


120 


-83 


25 


0.5 




pCO 


K/LS 


12.47 





240 


130 


47 


30 


0.25 




pCO 


K/LS 


14.4 




O 


-94 


51 


-114 


40 






pCO 


NTT/LS 


X 




64 


38 








21 


p 


K/N 


10.61 







583 


200 


327 


60 


0.1 




pCO 


NTT/LS 


11.73 




429 


105 


-105 


40 


0.1 


7SE 


p/early 


NTT/LS 


11.29 





3- 


-360 
-110 


80 
170 






0.5 
0.25 




p 
p 


Keck 
Keck 


12.35 
11.22 






a. 


-130 


40 


-230 


30 




34W 


p/early 


NTT/LS 


10.86 


2 


and 


115 


45 








3 


p/early 


NTT 


10.62 


2 






-90 


80 


0.25 


13SSE 


early 


LaSilla96 


X 













30 






early 


LaSilla96 


8.98 





-200 


80 






0.5 




p 


NTT 


12.68 


2 


O 






-15 


30 






CO 


LaSlll94 


12.89 


1 








-74 


30 


1 


13 W 


CO 


LaSIII94 


10.74 


2 








-98 


20 


1 




CO 


LaSIU94 


12.07 


1 








-103 


25 


1 




CO 


LaSIII94 


X 






225 


69 






0.25 




p 


NTT 


12.62 











-178 


50 


1 




CO 


LaSIII94 


10.68 











86 


25 






CO 


LaSIII94 












-300 


200 


0.25 


IS 


early 


LaSina96 


12.63 











17 


25 


0.5 


20 


CO 


LaSIU94 


10.52 







-118 


35 


-103 


15 


1 


7 


pCO 


NTT/LS 


6.37 


1 








32 


25 


0.5 




CO 


LaSm94 


12.66 


1 








-178 


40 


1 




CO 


LaSIU94 


9.49 


2 








170 


50 


1 


9NW 


early 


LaSilla96 


12.16 


2 








72 


30 


1 




CO 


LaSm94 


12.99 


1 








-40 


50 


1 


7E 


early 


LaSilla96 


X 










-68 


40 


1 




CO 


LaSIU94 


X 










107 


20 


1 


2S 


CO 


LaSIU94 


10.27 


2 
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CO 



> 

CO 

o 
o 
p 

o 
o 
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o 
o 
o 



Table 1 (cont'd 



R(SgrA*) 


Ax(SgrA*) 


Ay(SgrA*) 








AVy 


Vx AVx Vy AVy 


Vrf 


AVrf 


v„ 


AVjf 


Vj 


Av, 


weight 


name 


type 


source 


K-mag 


variabililty 


CD 


6.18 


(arsecs) 
5.70 


2.40 










all velocities are in km/s 










-34 


30 


1 


- 


CO 


LaSin94 


11.85 


2 


Itel 


635 


-3.91 


5.00 


75 


60 


-100 


80 




75 


60 


-100 


8« 


-300 


50 


1 


7W 


p/early 


NTT/LS 


11.85 


2 


fa 


6.62 


6.60 


0.50 




















29 


25 


1 


1NE<3) 


CO 


LaSin94 


10.86 


2 




6.72 


-6.20 


-2.60 




















-23 


20 


1 


- 


CO 


LaSm94 


12.41 


1 





6.77 


6.60 


1.50 




















-8 


25 


1 


1NE(2) 


CO 


LaSIll94 


X 




6.86 


-4.10 


-5.50 




















117 


40 


1 




CO 


LaSIll94 


11.69 


1 




7.34 


-2.50 


6.90 




















104 


40 


1 




CO 


LaSIll94 






P 


7.37 


7.30 


1.00 




















112 


25 


1 


INE(l) 


CO 


LaSIll94 


10.37 


1 


3 


7.47 


7.40 


-1.00 




















29 


25 


1 


ISE(l) 


CO 


LaSlll94 


10.73 





o 


7.52 


5.70 


-4.90 




















-79 


30 


1 


- 


CO 


LaSI1194 


X 




cc 


7.65 


-3.50 


-6.80 




















-96 


20 


1 


12N 


CO 


LaSIll94 


8.86 


2 


in 


7.80 


-6.30 


-4.60 


20 


60 


-45 


60 




20 


60 


-45 


60 


-97 


30 


1 




pCO 


NTT/LSS 


11.66 


2 


7.83 


7.00 


3.50 




















7 


30 


1 


lOEW 


CO 


LaSIll94 


X 




the G, 


7.94 


-0.86 


7.89 


-194 


90 


13 


48 




-194 


90 


13 


48 






1 


26 


p 


NTT 


11.06 





7.96 


-3.80 


7.00 




















-153 


80 


0.5 




CO 


LaSIll94 






8.03 


-7.10 


3.75 


-70 


80 


-120 


80 




-70 


80 


-120 


80 






0.5 




p 


NTT 


12.06 


2 


8.04 


-5.70 


5.67 


60 


60 


75 


30 




60 


60 


75 


30 


-77 


30 


1 


BHA4E 


pCO 


NTT/LS 


10.55 






8.10 


3.50 


7.30 




















-6 


50 


1 




CO 


LaSIll94 


11.51 


1 


act 


8.26 


-8.10 


1.60 




















-150 


70 


1 


6W 


early 


LaSilla94 


10.9 


2 


8.33 


-2.56 


7.93 


189 


85 


-52 


59 




189 


85 


-52 


59 






0.5 


86 


P 


NTT 


12.83 


1 


n 


8.34 


-7.20 


-4.20 




















150 


70 


1 


AFNW 


early 


LaSilla96 


X 




; ce 


8.42 


-1.90 


8.20 




















-33 


30 


1 


- 


CO 


LaSIll94 


X 




8.49 


2.20 


-8J0 




















88 


25 


1 


- 


CO 


LaSIll94 


12.19 


1 


C3 


8.59 


5.50 


6.60 




















77 


40 


1 


- 


CO 


LaSUI94 


X 






8.61 


8.40 


1.90 




















108 


30 


1 


- 


CO 


LaSDBM 


12.5 


2 


CD 


8.82 


0.60 


-8.80 




















19 


20 




14N 


CO 


LaSin94 


9.45 


1 




8.82 


-6.45 


6.02 


-270 


100 


175 


100 




-270 


100 


175 


m 


100 


30 


0.25 


BHA4W 


pCO 


NTT/LS 


10.44 


2 


1 


8.99 


6.00 


-6.70 




















-283 


25 






CO 


LaSDI94 


s 






9.04 


0.90 


9.00 




















23 


30 






CO 


LaSm94 


X 




-i 


9.06 


5.50 


-7J0 




















-300 


25 




9 


CO 


LaSin94 


8.94 


2 


9.09 


8.30 


3.70 




















-55 


30 




lOEE 


CO 


LaSm94 


11.56 


2 




9.45 


7.30 


-6.0« 




















-58 


35 






CO 


LaSIU94 


X 






9.53 


-0.70 


-9.50 




















29 


20 




14SW 


CO 


LaSIU94 


10.98 





o 


9.55 


-9.10 


2.90 




















-135 


40 






CO 


LaSm94 


X 




» su 


9.69 


-3.60 


-9.0fl 




















51 


25 


0.5 


12S 


CO 


LaSIU94 


9.94 


2 


9.73 


-9.15 


-3M 




















250 


100 




AFNWW 


early 


LaSilla96 


11.67 


2 




9.86 


-9.80 


-1.10 




















-50 


30 






CO 


LaSIU94 


X 






9.9fl 


9.30 


-3.40 




















22 


50 






CO 


LaSIU94 


X 




An 


9.91 


-7.30 


-6.70 




















140 


50 




AF 


early 


LaSilla94 


10.51 


2 


10.01 


3.70 


9.30 




















67 


40 






CO 


LaSm94 


X 




isotropy 



to 



Table 1 (cont'd 

1 2 3 4 5 6 7 8 9 10 II 12 13 

R(SgrA*) Ax(SgrA*) Ay(SgrA*) Av^ Vy Avy Av^ Vy Avy Vrf Av^ 

(arsecs) all velocities are in km/s 



@ 

o 
o 
o 
o 

> 

CO 



> 

CO 

o 
o 
o 

o 
o 
o 

i 

o 
o 



10.03 


-8.25 


-5.70 


10.06 


6.70 


-730 


10.08 


5.10 


8.70 


10.21 


-1.50 


10.10 


10.33 


2.60 


-10.00 


10.44 


-3.60 


980 


10.47 


-9.20 


-5.00 


10.51 


9.40 


-4.70 


10.60 


8.30 


-6.60 


10.61 


10.60 


0.50 


10.72 


5.50 


-9.20 


11.18 


7.60 


8.20 


11.24 


11.20 


-1.00 


11.30 


-5.80 


9.70 


11.32 


-8.20 


7.80 


11.38 


-8.10 


-8.00 


11.41 


-10.10 


-5.30 


11.47 


1.30 


11.40 


11.50 


-8.80 


7.40 


11.51 


1.60 


11.40 


11.52 


-2.70 


11.20 


11.61 


-11.00 


-3.70 


11.74 


9.50 


6.90 


11.83 


-0.90 


11.80 


11.87 


6.70 


9.80 


11.91 


4.80 


10.90 


11.98 


^.50 


11.10 


12.00 


0.20 


-12.00 


12.01 


11.90 


-1.60 


12.02 


11.50 


3.50 


12.05 


3.60 


11.50 


12J5 


11.80 


-3J0 


12.48 


5.50 


-IIJO 


12.59 


10.60 


-6.80 


12.66 


-9.20 


8.70 


12.71 


11.60 


5J0 


12.76 


3.40 


-12.30 


12.95 


12.50 


-3.40 


13.06 


-6.20 


11.50 


13.11 


6.10 


11.60 





















weight 


name 


type 


source 


K-mag 


170 


70 


0.5 




early 


LaSilla96 


11.9 


-283 


30 






CO 


LaSIII94 




147 


30 






CO 


LaSni94 




-230 


50 




15SW 




LaSilla94 




19 


30 


0.5 




CO 


LaSlll94 




134 


25 






CO 


LaSlll94 




-20 


70 




Blum 




LaSilla94 




82 


40 






CO 


LaSIll94 




-103 


35 






CO 


LaS11194 


12.77 


227 


25 






CO 


LaSIll94 




200 


50 




9S 




LaSilla94 




105 


35 






CO 


LaSlll94 




177 


30 






CO 


LaSlll94 




7 


30 






CO 


LaSlll94 




77 


30 






CO 


LaSlll94 




-18 


30 






CO 


LaSlll94 




87 


50 






CO 


LaSlll94 




17 


25 




15N 


CO 


LaSlll94 




59 


30 






CO 


LaSlll94 




-80 


50 




15NE 




LaSilla94 


11.4 


-68 


30 






CO 


LaSlll94 




-8 


40 






CO 


LaSIll94 




12 


50 






CO 


LaSIll94 


12.62 


45 


20 






CO 


LaSlll94 




37 


30 






CO 


LaSnijM 




121 


30 






CO 


LaSIII94 




-88 


50 






CO 


LaSni!M 




41 


30 






CO 






107 


30 






CO 


LaSIU94 


X 


137 


25 






CO 


LaSIU94 


s 


57 


30 






CO 


LaSIU94 


s 


27 


60 






CO 


LaSIU94 


X 


171 


35 






CO 


LaSIU94 


X 


-93 


25 




28 


CO 


LaSIU94 


X 


42 


50 






CO 


LaSIU94 


X 


114 


25 






CO 


LaSIU94 


X 


112 


20 




BHA17 


CO 


LaSIU94 


X 


347 


60 






CO 


LaSIU94 


X 


-23 


50 






CO 


LaSIU94 


X 


111 


30 






CO 


LaSIU94 


X 



S 

iS! 

O 



TO 

a. 

s 



@ 

o 
o 
o 
o 

> 

CO 



CO 

o 
o 
p 

o 
o 

? 

o 
o 
o 



Table 1 (cont'd 



p3 1 2 3 4 5 6 7 8 9 10 II 12 13 

> R(SgrA*) Ax(SgrA*) Ay(SgrA*) Av, Vy Avy v, Av^ Vy Avy Av^ 





(arsecs) 




13.15 


7.50 


-10.80 


13J2 


13J0 


-0.70 


13.28 


8.50 


10.20 


13.30 


9.10 


-9.70 


13.34 


11.70 


6.40 


13.49 


13.20 


2.80 


13.50 


-0.20 


-13.50 


13.66 


-1.30 


13.60 


13.67 


10.20 


9.10 


13.87 


5.10 


-12.90 


13.88 


-13.10 


-4.60 


13.90 


11.50 


-7.80 


13.98 


12.90 


-5.40 


14.04 


6.40 


-12.50 


14.06 


13.10 


5.10 


14.06 


-7.80 


-11.70 


14.06 


-11.70 


7.80 


14.25 


14.20 


-1.20 


14.37 


-5.90 


-13.10 


14.38 


-11.60 


-8.50 


14.48 


3.70 


14.00 


14.58 


-10.80 


9.80 


14.64 


-8.80 


-11.70 


14.64 


5.90 


13.40 


14.65 


10.70 


10.00 


14.71 


-14.00 


4.50 


14.88 


-5.80 


13.70 


15.03 


2.00 


-14.90 


15.08 


6.30 


-13.70 


15.11 


13.00 


-7.70 


15.19 


11.60 


-9.80 


15ai 


-3.50 


-14.80 


15 J5 


-13.80 


6.50 


15J1 


15.30 


0.40 


1535 


-7.30 


13.50 


1537 


^.50 


-14.70 


15.59 


-9.70 


12.20 


15.78 


14.90 


5J0 


15.81 


-13J0 


8.70 


15.8« 


-10.60 


-11.80 



all velocities are in Icm/s 



16 


17 




19 


20 


21 


Vz 


Av, 


weight 


name 


type 


source 


18 


40 






CO 


LaSIU94 


187 


30 






CO 


LaSII194 


97 


30 






CO 


LaSIU94 


-138 


70 






CO 


LaSIll94 


67 


40 






CO 


LaSIll94 


106 


40 


0.5 




CO 


LaSIll94 


-38 


50 






CO 


LaSIll94 


-45 


50 






CO 


LaSIll94 


-56 


40 






CO 


LaSI1194 


137 


40 






CO 


LaSIll94 


-93 


50 






CO 


LaSIll94 


-163 


30 






CO 


LaSIll94 


-36 


50 






CO 


LaSIll94 


87 


35 






CO 


LaSIll94 


185 


25 




17 


CO 


LaSIll94 


-73 


40 






CO 


LaSIll94 


130 


60 






CO 


LaSIll94 


-183 


50 






CO 


LaSIll94 


-58 


70 






CO 


LaSIll94 


72 


25 






CO 


LaSIll94 


117 


60 






CO 


LaSIll94 


77 


45 






CO 


LaSIll94 


-45 


40 




BHA2 


CO 


LaSIll94 


-83 


70 






CO 


LaSIll94 


67 


40 






CO 


LaSnilM 


33 


50 






CO 


LaSII194 


-48 


40 






CO 


LaSIU94 


64 


25 






CO 




31 


50 






CO 


LaSIU94 


-38 


40 






CO 


LaSm94 


34 


35 






CO 


LaSm94 


67 


60 






CO 


LaSIU94 


32 


50 






CO 


LaSIU94 


-73 


70 






CO 


LaSm94 


27 


40 






CO 


LaSIU94 


-33 


70 






CO 


LaSIU94 


-24 


30 




IISW 


CO 


LaSm94 


121 


30 






CO 


LaSIU94 


3 


35 






CO 


LaSIU94 


-63 


40 






CO 


LaSIU94 



K-mag variabililty 



(1) 

D 

3 

o 

cc 



cr 

(1) 

Q 

o 



O 

(1) 

I 



o 
o 

Co 



Co 
O 

I 



Table 1 (cont'd 

1 2 3 4 5 6 7 8 9 10 II 12 13 

R(SgrA*) Ax(SgrA*) Ay(SgrA*) Av^ Vy Avy Av^ Vy Avy Vrf Av^ 

all velocities are in km/s 



@ 

o 
o 
o 
o 

> 

CO 



> 

CO 

o 
o 
o 

o 
o 
o 

i 

o 
o 





(arsecs) 




16.01 


-0.60 


-16.00 


16.04 


2.10 


15.90 


16.10 


6.10 


14.90 


16.18 


2.40 


-16.00 


16.18 


12.80 


-9.90 


16.20 


11.80 


-11.10 


16.30 


4.40 


15.70 


16.35 


-8.60 


13.90 


16.35 


4.20 


-15.80 


16.46 


-16.40 


1.40 


16.46 


14.50 


-7.80 


16.49 


10.40 


-12.80 


16.65 


-16.00 


460 


16.77 


-13.00 


-10.60 


16.95 


-3.90 


-16.50 


16.98 


-9.00 


-14.40 


17.01 


-16.90 


1.90 


17.42 


-6.40 


16.20 


17.45 


1.30 


-17.40 


17.49 


-15.50 


8.10 


17.54 


8.40 


-15.40 


17.66 


-1.50 


17.60 


17.98 


-6.40 


-16.80 


18.02 


-8.10 


16.10 


18.08 


-7.40 


-16.50 


18.10 


0.20 


-18.10 


18.15 


-9.60 


-15.40 


1800 


-9.20 


15.70 


18J8 


11.00 


-14.60 


18J8 


9.70 


-15.50 


18J0 


0.20 


-18.30 


18J4 


-5.80 


17.40 


18.51 


0.50 


18.50 


18.54 


-3.00 


-18J0 


18.58 


15.60 


-10.10 


18.66 


-17.00 


7.70 


18.81 


13.40 


-13.20 


18.92 


17.50 


-7J0 


19.01 


-5.50 


-18J0 


19.03 


-10.60 


-15.80 



ight 



type 



source K-mag variabililty 



S 

iS! 



19 


40 






CO 


LaSin94 




57 


50 






CO 


LaSm94 




92 


50 






CO 


LaSm94 


^ . 


-28 


60 


1 




CO 


LaSm94 


hon, . 


50 


40 






CO 


LaSIll94 


81 


35 






CO 


LaSIll94 


-36 


60 






CO 


LaSIll94 




14 


30 




IINE 


CO 


LaSI1194 


.Ecka 


-151 


80 






CO 


LaSI1194 


-60 


60 


0.5 




CO 


LaSIll94 


50 


60 






CO 


LaSIll94 




-138 


50 






CO 


LaSIll94 


-129 


50 






CO 


LaSIll94 




-103 


40 






CO 


LaSIll94 


TO 


-148 


50 






CO 


LaSIll94 


37 


40 






CO 


LaSIll94 


rha 


97 


40 






CO 


LaSIll94 


-118 


70 






CO 


LaSIll94 


rd an 


70 


50 






CO 


LaSIll94 


163 


25 






CO 


LaSIll94 


27 


60 


0.5 




CO 


LaSI1194 




107 


30 






CO 


LaSI1194 


T. Ot 


-213 


60 






CO 


LaSI1194 


-29 


50 


0.5 




CO 


LaSIll94 


-84 


50 


0.5 




CO 


LaSIU94 




17 


30 






CO 


LaSm94 




-63 


40 






CO 


LaSm94 




-8 


50 


0.5 




CO 


LaSm94 




17 


70 


0.5 




CO 


LaSIU94 




107 


70 


0.5 




CO 


LaSm94 







30 






CO 


LaSm94 




-8 


60 






CO 


LaSIU94 




-53 


60 






CO 


LaSIU94 




-203 


50 




BHA7 


CO 


LaSm94 




2 


40 






CO 


LaSIU94 




87 


40 






CO 


LaSIU94 




7 


50 






CO 


LaSm94 




-93 


60 






CO 


LaSIU94 




-63 


80 






CO 


LaSIU94 




157 


70 






CO 


LaSIU94 
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Table 2 

Velocity Dispersions, Anisotropies and Projected Mass Estimators 



Annulus 


N 


<R> 




Or 




P 


M(BT) 


M(VT) 


M(LM) 




(stars) 


(arcsecs) 


(km/s) 


(kiWs) 








(10' M^) 




0<R<0.8" 


17 


0.3 


334(58) 


503(88) 


0.66(0.15) 


0.66(0.2) 








0<R<1" 


19 


0.4 


347(57) 


479(79) 


0.72(0.14) 


0.58(0.23) 


2.2(0.5) 


2.2(0.7) 


2.4(0.8) 


1<R<2" 


26 


1.4 


243(34) 


168(24) 


1.45(0.29) 


-3.6(6.1) 


2.8(0.5) 


2.7(0.8) 


2.5(0.7) 


2<R<3" 


26 


2.3 


175(25) 


135(19) 


1.30(0.26) 


-1.54(2.3) 


3.2(0.6) 


2.8(0.8) 


2.8(0.8) 


3<R<5" 


23 


3.5 


153(22) 


148(22) 


1.03(0.21) 


-0.11(0.72) 


3.6(0.7) 


3.0(0.9) 


3.6(1.1) 


5<R<9" 


9 


7.3 


115(28) 


117(28) 


0.98(0.34) 


0.06(0.9) 


3.7(1.2) 


3.8(1.8) 


3.6(1.7) 


0<R<5" 


95 


1.4 


236(17) 


250(18) 


0.94(0.1) 


0.15(0.23) 


3.06(0.31) 


2.62(0.38) 


2.91(0.42) 








243(18) 


246(18) 


0.99(0.1) 


0.04(0.29) 


3.14(0.32) 


2.46(0.36) 


2.93(0.42) 



• Column 21 lists the source(s) of measurement for the specific data point (Keck (K) and NTT (N) for proper motions, 
LaSilla (LS), or 'aU'). 

• Column 22 lists the K-magnitude of the star, and column 23 makes a statement on its variability. If available, we used 
the K- magnitudes from the comprehensive variability study of Ott et al. (1999), otherwise we list the magnitudes of Chez et 
al. (1998), corrected for 5mA'=-0.4 to account for a small calibration offset between the Ott et al. and Chez et al. sets. If the 
value in the last column is 0, the star was not or not significantly variable in the 1992 to 1998 monitoring campaign of Ott et 
al. A value of 1 indicates that the star showed a statistically significant but weak variability. A value of 2 indicates that the 
star was strongly variable in the Ott et al. observations. 



3 KINEMATICS OF THE GALACTIC CENTRE STAR CLUSTER 

The velocity determinations in Table 1 are significantly improved over our earlier work and over Ghez et al. (1998). Many light- 
of-sight and sky velocities now have errors less than 50 kms~^, with the best velocity determinations (±20 to 25 km/s) mainly 
limited by systematic effects (e.g. establishing reference frame from the moving stars themselves and removing distortions 
in the imaging), rather than by statistical errors (positional accuracy of stellar positions). Of the 2 x 104 proper motions in 
Table 1, 48 (23%) are determined to 4(t or better. 5 proper motions are determined at the >10a level. Of the 227 line-of-sight 
velocities 38 (17%) are determined to or better. For 14 (of 29) Hel emission line stars and for 18 late type stars we now have 
determinations of all three velocity components. With this improved data set it is now possible to investigate in more detail 
the kinematic parameters of individual stars and/or small groups of stars. In order to remove as well as possible measurement 
and calibration bias and zero point offsets, we subtract in all of our calculations below for each velocity measurement i the 
mean velocity of the sample (v) linearly and the velocity uncertainty error (vi) in squares when computing velocity dispersions 
etc., 

(«?)intrmsic = ({"i - (tl) }^ )mcasm-od " CrrOr^ (Vi) . (1) 

3.1 Tests for Anisotropy 

As proposed by Eckart and Genzel (1996, 1997), a first simple (but coarse) test for anisotropy in the data (and/or a Sun- 
Galactic centre distance significantly different from Rq = 8kpc) is to compare the sky and line-of-sight velocities of individual 
stars. Fig. (^ is a plot of 7pz = (vp^op - v^)/ (vp^^p -I- v^) as a function of projected separation R for the 32 stars with three 
measured velocities. Here Vprop is the root mean square of the x- and y-sky motions, Vz is the line-of-sight motion. In this plot 
stars with = —1 have Vz 3> fprop stars with 7p2=-|-l have v^, Vprop- Fig. (^) shows no obvious sign for such an anisotropy. 
This is probably not surprising as the line-of-sight and sky velocities both contain linear combinations of the intrinsic radial 
and tangential components of the velocity ellipsoid. The result in Fig. (^ that the sample expectation value for the proper 
motion velocity dispersion is the same (within statistical uncertainties) as the line-of-sight dispersion, < Vp^op >—< >, is 
consistent with the assumption that we are observing a spherically symmetric cluster. The virial theorem guarantees that this 
results holds independent of internal anisotropy (see equations (^) (^) (^)) . 
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Figure 4. The anisotropy parameter = (Vpj.gp — v^) / (Vp^op + v^) for all 32 stars with all three velocity components measured. Error 
bars are determined through error propagation from the velocity uncertainties in Table 1. 



To investigate intrinsie kinematic anisotropics it is therefore necessary to explicitly decompose the observed motions into 
projections of the intrinsic velocity components. Assuming that the velocity ellipsoid of a selected (sub-) sample of stars 
separates in spherical coordinates and denoting the components of velocity dispersion parallel and perpendicular/tangential 
to the radius vector r as ar and at, the line-of-sight component az{R,z) is then given by 

= crr{r) cos^ 6 + <jl{r) sin^ 9, (2) 

where cos^ = r • z/r and z is the unit vector along the line-of-sight. The components of the velocity dispersion parallel (R) 
and perpendicular (T) to the projected radius vector on the sky R are given by 

= cr^(r) sin^ -h CTj (r) cos^ 0, (3) 

and 

4 = a?(r) . (4) 

Given the spatial density distribution n(r) of the selected sample of stars (assumed to be spherically symmetric) the line- 
of-sight averaged, density weighted value of the projected radial velocity dispersion of the sample at R, an, can then be 
computed from the relationship 

m)crUR) = P n{z){al{r) sin^ 6 + erf (r) cos^ e)Az = 2 [al{r)(R/rf + a?(r)(l - [R/rf] ^J,'^^^'^y,, (5) 
where S(_R) is the stellar surface density at R, 

-(«)-/:^ 

Similar equations hold for a^{R) and crj(i^). The global expectation value of the projected radial dispersion is given by 

(aji) = {2n/N) / n(r)r^((T^(r) sin^ 61 (T?(r) cos^ 61) sin6>d6idr = 2/3(ct^) -|- l/3(cTf) , (7) 

Jo Jo 

where N is the number of stars in the selected sample. Likewise one finds 

(al) = {2-k/N) / / n{r)r^{al{r) cos^ 9 + a? (r) sin^ 6*) sin6id6'dr = l/3(c7^) + 2/2,{al), (8) 
Jo Jo 



and 



Wt) = {c^h (9) 
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Figure 5. The anisotropy parameter 7Tfl = (fj? — ^^)/('''7' + ^ji) for the 12 Hel stars within R < 5" (top left), for the 14 late type 
stars within i? < 5" (bottom left), for all stars inside 5" with proper motions (bottom right) and for stars within R < 0.8" of SgrA* (top 
right). Error bars are determined through error propagation from the velocity uncertainties in Table 1. 



(Leonard & Merritt 1989). Deviations of the velocity ellipsoid from isotropy are commonly expressed in terms of the anisotropy 
parameter f3 = 1 — at /a^- Its globally averaged value is given by 

(13) = 1- {a^,)/{al) = 3((a^) - {a^)/(^al) - (4» • (10) 

An isotropic cluster ((/3)=0) has (ar) = (crt) and {(JR} = {aT) ■ A cluster with only radial orbits ({/3) = 1) has {at}={o'T)=0 or 
(<7"i{) » (ax)- A cluster with only tangential orbits {((3) = -oo) has (ar) =0 and (crR)= (crT)/\/3. Thus radial anisotropy is 
easier to see in the proper motions than tangential anisotropy. Table 2 gives the values of (ctr), (ctt), {crR)/(aT) and (/?), 
computed for all stars with proper motions and for different ranges of projected radii from SgrA*. Errors in those quantities 
are derived from statistics and error propagation. Below we use Monte Carlo simulations to investigate the uncertainties in 
the derived anisotropy parameters more throroughly. The proper motions of the entire sample of stars as well as the stars in 
the range R>3" are consistent with isotropy. At R=l to 3" there is a (marginal) trend for the stars to be more on tangential 
orbits. In the central arcsecond the stars on average appear to be on radial orbits. The statistical significance of this departure 
from isotropy for the 17 stars at R<0.8" from SgrA* appears to be 3.3cr in terms of propagated errors for j3 (Table 2); however, 
the Monte Carlo simulations of Section 4.3 below show that the distribution of (3 is very broad and the isotropic /? = is still 
within somewhat more than Icr equivalent for such a small sample. Excluding the faint stars in the SgrA* cluster at R< 1" 
the remaining proper motions between R=l" and 5" deviate from isotropy in the direction of tangential orbits at the 2a level. 
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Figure 6. Possible orbits for SI and S2, the stars closest to SgrA*. The positions of Sl(left) and S2 (right) between 1992 and 1999 and 
their uncertainties are shown as crosses. Three possible bound orbits in the potential of a 2 to 3 X 10^ Mq point mass are shown. Orbit 
A (continuous curve) has the largest possible separation from SgrA* and orbit B (short dashes) the largest line of sight velocity. The 
assumption that SI and S2 are in the same plane of the sky as the central mass and have no line of sight velocity results in the orbit 
with the largest curvature (orbit C, dash-dotted). 



A second and more sensitive test for anisotropy is a comparison of the projected radial iyB.) and projected tangential 
(vt) components of the sky velocities of individual stars. Fig. (1) gives plots of 7th ~ («t ~ ^r)/{''^t + ^r) for different 
selections of our data. Considering all 104 proper motions, Fig. (E) (bottom right) indicates a fairly even distribution of jtrs, 
without obvious overall bias indicating anisotropy, perhaps a slight predominance of tangential orbits (compare Fig. The 
same is true if only the late type stars of the proper motion sample are considered (Fig. ( ^) , bottom left) . 

A different and fairly clearcut picture emerges when one considers the (much younger) early type stars. Fig. (1) (top left) 
clearly indicates that with one exception all bright Hel emission line stars within R=5" are on projected tangential orbits 
(7Tfl~+l) and therefore (see the discussion after equation (^o|) ), largely on true tangential or circular orbits. In contrast, 
more than half of the faint (mK«13 to 16) stars within 1" of SgrA* (SgrA* cluster) are predominantly on radial {"/tr < 
0) orbits (top right panel of Fig. (^)). We conclude that the early type stars m our proper motion set do show significant 
anisotropy.. We will show below that the main cause of the tangential anisotropy is a global rotation of the early type stars. 



3.2 Orbits for the innermost stars 

We have also modeled the orbits of several of the individual fast moving stars in the SgrA* cluster. As an example we plot in 
Fig. (^ the measured 1992-1999 NTT positions of SI and S2 with respect to SgrA*, along with the projection of a few possible 
trajectories. The three plotted orbits represent extreme choices of the orbital parameters in the potential of a central compact 
mass. For orbit A we assumed the largest possible current separation from SgrA* for bound orbits with Vz = 0. For orbit B 
we took the largest line of sight velocity at z = under the boundary condition that SI and S2 are still bound to SgrA*. The 
assumption that S1/S2 are in the same plane of the sky as the central mass and have no line of sight velocity results in the 
orbit with the largest curvature (orbit C). Although no unique orbit can yet be determined from the data, our analysis shows 
that most of the high velocity stars in the SgrA* cluster can be bound to a central mass of 3 x W^Mq with a distribution 
of line of sight positions and velocities that is consistent with the projected dimensions and velocity dispersion of the cluster. 
Orbits with radii of curvature comparable to their projected radii from SgrA* (orbit C) can already be excluded for these stars. 
SI, S2 and several other fast moving stars around SgrA* must be on plunging (radial) orbits or on very elliptical/parabolic 
orbits with semi-axes much greater than the current projected separations from SgrA*, as already indicated by the analysis of 
the velocities in the preceding paragraph. The stellar orbits in the central cluster will not be simple closed ellipses. Especially 
for orbits with large eccentricities, node-rotation due to the non-Keplerian potential in the extended stellar cluster as well 
as relativistic periastron rotations will make the trajectories for the individual stars 'rosette'-like. A more detailed analysis 
of orbits has to await longer time baselines for the proper motion measurements and the detection of orbit curvature {— 
acceleration), as well as measurements of the line of sight velocities of the central stars. 
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3.3 Anisotropy and Relaxation Time 

These deviations from isotropy for the early-type stars are consistent with their young ages as compared to the relaxation 
time. Within the central stellar core, the two-body relaxation time for a star of mass mio (in units of 10 Mq) is given by 

trim) = 10^''V?oo/(p6.6mio{?niV./13}) (yrs). (11) 

Here pe.e is the density of the nuclear star cluster (in units of 4 x IO^Mqpc"^) and crioo is the velocity dispersion in units 
of 100 km/s within the core radius Qof 0.38 (-H0. 25,-0. 15) pc (Genzel et al. 1996). N. is the number of stars in the core 
(~ 3to5 X lO""). The life-time of the upper main sequence phase scales approximately as tms ~ lO^'^m^Q^'^yr and the duration 
of the red-/blue-giant or supergiant phases is typically 10 to 30% of tms (e.g. Meynet et al. 1994). The ratio of relaxation 
time to stellar life time thus is 

X = tr/tms~2Am°J. (12) 

A number of authors have shown have shown that the Hel emission line stars are high mass (30 to 120 Mq), post main 
sequence blue supergiants (Allen et al. 1990, Krabbe et al. 1991, 1995, Najarro et al. 1994, 1997, Libonate et al. 1995, Blum, 
Sellgren and dePoy 1995 a,b, Tamblyn et al. 1996, Ott et al. 1999). Their ages range between t « 2and9 x 10® yrs (Najarro 
et al. 1994, 1997, Krabbe et al. 1995). The massive stars are probably the 'tip of the iceberg' of a component of young stars 
of total mass ^10" Mq that were formed a few million years ago in an extended starburst episode (Krabbe et al. 1995). 
The massive stars are somewhat older than their main sequence age. Their main sequence life time is much greater than the 
dynamical time scale (tms 3> tdyn — lO'^yr) but they have not had time to dynamically relax through multiple interactions 
with other stars [x^ !)• Their present kinematic properties thus reflect the initial conditions with which they were born and 
the starburst must have been triggered near their present orbits. The situation is different for the observed late type stars. 
For M-giants of mass 1.5 to 3 Mq (and ages >10^ years) and for luminous asymptotic giant branch (AGB) stars of mass 3 
to 8 Mq (and ages >10* years), x is comparable to or smaller than unity. Such stars should have had sufficient time to be 
scattered and relax in the central potential. 



4 MONTE CARLO SIMULATIONS 

Because the velocity measurement errors are often large and the number of measured velocities is still relatively small, we 
need to investigate the expected errors in the velocity anisotropy in more detail to get a more quantitative estimate whether 
the observed anisotropy is statistically significant. Therefore we now describe theoretical 'measurements' on Monte Carlo star 
clusters with comparable numbers of stars. The next subsection describes the models from which the artificial data are drawn. 



4.1 Anisotropic distribution functions 

We construct some simple anisotropic, scale-free spherical distribution functions f(e,h) for stars with specific energy e and 
specific angular momentum ft in a potential ip. These are computed from the formula (for a derivation and its generalization 
to non integer index n, see, e.g., Pichon & Gerhard, in preparation): 

f[e,h) = (-e)"+^__^2^ / d y+^^,„+4^^^ 



47rV^r(| +n) \Ar'^ip 



(13) 

where specific energy and angular momentum are given by 

e — ^{v^. + vl) — ^ and h — rvt- (14) 

Neglecting the self gravity of the star cluster in the vicinity of the central black hole we write 

GM f r \ -5/2 

'(/>(r) = , and p(r) = po — , (15) 

r Vro/ 

where M is the black hole mass and G the Gravitational constant. The slope of the cluster density proflle is chosen to provide 
a compromise between the observed m(K)<15 number counts (see below) and the observed distribution of the innermost 
SgrA* cluster stars. 

The resulting distribution from Eqs. (p^)-([l5|) reads 

r, ,^ , r+^/,2n-i _ 3^ 35 (-e)^^ 231(-£)j/i' 1287 (-e)^ ft' 46189 {-e)i h' 96577 (-e)^ ft' . . 

f[e,n) cc [ e) h 64ft7r3' 1287r3 ' 2567r3 ' 5127r3 ' 71687r3 ' 61447r3 ^ ' 

oc {-e)[h/h4e)f"-^ (17) 



The core radius is defined here as the radius where the stellar density has fallen to half its central value. 
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Figure 7. left panel: Sky projection of 250 proper motions drawn from a star cluster with p oc r~ and /3* = 1/2 (radial anisotropy) 
and right panel: /3* = —1/2 (tangential anisotropy) in the potential of a point mass. The distinction between radially and tangentially 
anisotropic clusters is clearly visible. Radially anisotropic clusters are more easily recognizable on the basis of the proper motion patterns 
than tangentially anisotropic clusters. Note that it may appear incorrectly that most stars are moving outwards because the inner 
plunging orbits have longer vectors and seem to be moving radially outwards as the arrows overshoot the centre. 



for n = 0, 1, • • • 5. Here hc{e) denotes the angular momentum of the circular orbit at energy e. The units are such that the 
total mass and scale-radius, ro, of the star cluster are unity and GM = 1 . For this simple scale-free cluster in a Keplerian 
potential the distribution functions are also derived in Sections 2.2 and 3.1 of Gerhard (1991). 
The velocity dispersion corresponding to Eq. (jl^) is 



d 



I 2 2n+6\ 



2n+4 

r /9, 



and together with the Jeans equation 



[par] + 



d_ 

At 

this implies 



2/3 2 
— po-r = 
r 



pGM 



= 1/2 ■ 



(18) 



(19) 



(20) 



where the * subscript refers to the fact that this is the intrinsic anisotropy of the model. These simple models are therefore 
scale- free and have constant anisotropy parameters /?*, which in the following will be chosen to match the range of values 
found for the Galactic centre data set. 



4.2 Monte Carlo star clusters 

We generate A*' stars sampled regularly in radius with a cumulative mass profile corresponding to Eq. (^^. These are also 
required to obey Eq. (p^, i.e. the number of stars at radius r within dr with radial velocity Vr within dvr and tangential 
velocity vt within dvt is given by 

dA'' = 8n^ r^Vt f{\(v1 + v1) — ^^{r) ,rvt)Ardvtdvr . (21) 

Given a triplet (r, Vr,vt), we generate a random position vector r = r{cos[6) cos(<jf>), cos(S) sin((^), sin(S)} where z is along the 
line-of-sight as before, (j) is a. random number uniform in [0, 27r[ and sin(^) is a random number uniform in [—1, 1[. We also 
construct vg — Vt cos(x) and = vt sin{x) where X is a random number uniform in [0, 27r[. The velocity vector then reads 
V = VrBr + vgBe + v^G^. It is then straightforward to project the components of r and v onto the plane of the sky. 

Fig. ^ displays sky projections of the proper motion vectors of 250 stars drawn from = 1/2 and —1/2 clusters, 
respectively. The length of each arrow is proportional to the magnitude of the projected proper motion. The figure shows that 
a radially anisotropic cluster is readily recognized by the many stars with radial proper motions: intrinsically radial orbits 
remain radial when projected onto the sky, and the number of radial proper motions is a good indicator of the number of 
radial orbits. By contrast, intrinsically tangential orbits may appear tangential or radial in the sky plane, depending on the 
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Figure 8. Histograms of 7Ti{ = (i)^ — v'^)/(v^ "'"''I;) ^'-'^ ^ total number of 25 stars drawn from tlie simulations described in Section ^ 
Dark histograms show the mean number of stars per bin averaged over 50 draws of 25 stars each out of a 5000 stars cluster. The added 
light histogram shows the mean relative errors of the bin values. From top to bottom and left to right: /3t = 0.5, 0, —0.5, —3.5. Note that 
the isotropic model shows more stars on projected tangential than on projected radial orbits. 




Y 



Figure 9. Left panel: Theoretical ■ytr probability distribution (pdf) for star clusters with constant anisotropy given by /3* = 
1/2,-1/2,-3/2,-5/2,-7/2, —9/2, corresponding to expectation values for 7tr ofn/(n + l) =0, 1/2, 2/3, 3/4, 4/5, and 5/6 and variance 
of (1 + 2n)/(l + n)^ /(2 + n) = 1/2, 1/4, 5/36,7/80,3/50 and 11/252. These distributions are very non Gaussian and peaked near ±1. 
The /3* = 1/2 curve is symmetric (rather than that for the isotropic model) because we have defined ^tr in terms of the total on 
the sphere rather than one half that. For these pdfs the mean and standard deviation are well— defined for all values of /3*. Right panel: 
Projected ■jtr pdf for 2000 stars in clusters with /3* = 1/2,0, —1/2, —7/2 (dashed, full, dotted, and long dash-dotted lines). As for ■ytr, 
none of these curves is uniform. Note that all pdfs have significant tails for -Jtr = —1 because even purely tangential orbit may project 
to projected radial orbits on the sky. Most of the difference between the various models is in the relative number of stars in the radial 
and tangential peaks near ±1: The most tangentially anisotropic model (long dash-dotted line) has the most stars for Jtr G [0-5, 1[, 
whereas the radially anisotropic cluster (dashed line) is overabundant in near radial proper motions. The curve for the isotropic cluster 
is nearly symmetric. 



orientation of their orbital planes. Correspondingly, the projection of the tangentially anisotropic cluster in Fig. ([|), while 
showing fewer radial and more tangential proper motion vectors, still contains a significant number of the former. Therefore, 
tangentially anisotropic clusters are more difficult to recognize and discriminate from each other in terms of their apparent 
proper motion distributions. Once the model is sufficiently tangential, the ratio of radial to tangential proper motions is 
largely determined by the projection rather than by the intrinsic anisotropy. 
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Figure 10. Median (tliick line) and first and third quartiles (tiiin lines) for the distribution of (/3) values determined by Eq. (|lo[), from 
simulated proper motion samples, as a function of sample size (per bin) and for several values of true /3* for the underlying star cluster. 
The asymmetric width of the confidence bands reflects the asymmetric nature of the anisotropy parameter f}. For small samples Eq. ( [lo[ ) 
allows unphysical derived values for /3. Note the drift of the median as a function of which also reflects the fact that the pdf of j3 has 
inherited that of j3, with both skewed towards negative values. See Eq. ( |23| ) and Fig. The rather slow convergence of the estimator 

as a function of the number of stars for models with more negative f) is striking; this follows indirectly from Eq. (^sj) which has no 
well-defined moments. 



4.3 Anisotropy estimators 

From the Monte Carlo sample of stars we can estimate the previously used anisotropy indicators 7Ti?., and (Eq. ([ic| ) ). The 
histograms of the estimated 7Tfl in Fig. (^ confirm the above discussion quantitatively: (i) They show that radially anisotropic 
models are more easily recognized by their proper motion anisotropy than tangentially anisotropic models. Nonetheless, 
strongly tangentially anisotropic clusters are recognizable in terms of their many stars with "fxR near +1. (ii) The distribution 
of 7Tfi is slightly skewed towards positive values even for near- isotropic clusters, (iii) The histograms are always bimodal, i.e., 
have peaks near 7tjj = ±1. This is also recognizable in the data; cf. the bottom right panel of Fig. (^. 

These histograms are discrete realizations of the probability distribution for 7Tfl, and this in turn derives from the 
marginal probability distribution for the intrinsic quantity 'ytr ~ {'Vt — 'Vr)/{vr + vl), i.e., the number of stars with ^tr in a 
small interval A'ytr- Once the distribution function is known this is straightforward to compute: 
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Figure 11. Left Panel: Theoretical (3 probability distribution (pdf) for = 1/2, —1/2 ■ ■ ■ — 9/2 as labelled. These distributions are very 
skewed, and their means and standard deviations are ill-defined. Right Panel: Monte Carlo estimate of the probability distribution of $ 
for models corresponding to /3* = 1/2, 0, —1/2 and —7/2. The distributions result from applying Eq. ( |lo| ) to the values of (f^) and 
obtained for many random realizations of 20 stars from the respective star cluster models. 



which yields (after normalization) 



pdf(7t.)d7t, = '^5^^ d7t. where (n)!! = n(n - 2)(n - 4) • ■ ■ 1 . (22) 

TT (2n - l)!!vl - 7tr 

This pdf is illustrated in Fig. (^; it is strongly non-Gaussian and skewed for both n > and for an isotropic cluster. The 
reason why the isotropic curve is not symmetric is because we have defined ■jtr in terms of the total on the sphere rather 
than one half that. The main point of this diagram is the non-uniform and sometimes bimodal shape of the distribution. The 
distribution of the observed jtr pdf after projection is also shown in in Fig. (^ as a histogram for 5000 stars. Relative to 
the intrinsic pdf the number of (projected) radial orbits has been boosted, as discussed above; the distribution is now always 
bimodal. Thus in diagrams like Fig. ( |^ we should expect to always find an overabundance of stars near -yxR = ±1 compared 
to values near -yrR — 0, with the ratio N{'yTR ^ '^)/N{jtr ^ —1) containing the information about anisotropy. 

Fig. (^^ shows the median and first and third quartiles for the distribution of {{3) values determined by Eq. ( p^ from 
simulated proper motion samples, as a function of sample size and for several values of true /?* for the underlying star cluster. 
These confidence bands are especially wide for negative values of /3* because it is a very asymmetric indicator of anisotropy. 
Indeed the marginal propability distribution for the 13 values as determined from individual stellar velocities is given by 
(following the derivation of Eq. (p2|)) 

pdf{f3)df3 = ^^^!^I%^IHZfl^d/3 . (23) 
^■'^'"^ 7r(2n - 1)!!(3 - 2/3)"+i ' ^ ' 

This pdf is illustrated on the left panel of Fig. (pT[). Eq. ( ^) does not have any moments (i.e. the pdf does not fall off 
fast enough as a function of /? to allow for, say the mean and the variance to be computed). This implies that any estimator for 
its central value will be unreliable. The pdf for estimated via Monte Carlo simulations, has inherited these asymmetries; 
see the right panel of Fig. ( |l^) . Because of the observed skewness we expect the mean and the median to overestimate the 
anisotropy, especially for more negative models, as was indeed seen in Fig. (p^. 

We are now in a position to discuss the inferred anisotropies of the Galactic Centre star cluster (Fig. (|^ and Table 2) 
in more detail. Comparing with the distributions in Fig. (|^), the evidence for radial anisotropy in the central 0."8 rests on 
the absence of stars with jtr — 1, and the case for the tangential anisotropy of the Hel stars on the absence of stars with 
7th — —1. Based on the Monte Carlo models the evidence for anisotropy of the orbits is fairly solid. With larger proper motion 
samples, it may be best to compare with these distributions directly to estimate the anisotropy. The values for /3 estimated 
from Eq. (|l^, on the other hand, are quite uncertain. With this estimator being a quotient of observable dispersions, its 
distributions are very broad (Fig. (0)). The Monte Carlo simulations indicate that a sample of 500 stars with v/a(v) greater 
than 3 will be required in order to determine even /?* = 1/2 to an accuracy of ±0.2. 
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Figure 12. Line of siglit velocities of all Hel stars in Table 1, as a function of Dec-offset from SgrA*. Error bars a ztlu. 



In summary, the number of observed stars and the quahty of the derived velocities is already sufficient to state with 
some certainty that anisotropies in the orbits of (early type) stars are indeed present. To be consistent with the observed 
distribution of 7Tfl, the model clusters (assuming sphericity and cylindrical symmetry of the velocity ellipsoid) require fairly 
strong radial anisotropy at small radii, and tangential anisotropy for larger radii. However, the data are not yet suited to 
place accurate quantitative constraints on the anisotropy parameter (3 and its radial dependence. 



5 GLOBAL ROTATION OF THE EARLY TYPE CLUSTER 

As a group, the early type stars {= the starburst component) exhibit a well-defined overall angular momentum. The line- 
of-sight velocities of the 29 emission line stars follow a rotation pattern: blue-shifted radial velocities north, and red-shifted 
velocities south of the dynamic centre ( Fig. (^)). The apparent rotation axis is approximately east-west, within ±20°. The 
early type stars thus are in a counter-rotation with respect to general Galactic rotation, the latter showing blue-shifted material 
south and red-shifted material north of the Galactic centre. The rotation is fast (average «150 km/s) and is consistent with 
a Keplerian boundary for a 2 to 3 million solar mass central mass (Fig. ([l^)). Our results confirm and improve the earlier 
conclusions of Genzel et al. (1996). Note that the late type stars also show an overall rotation, but that is slow (few tens of 
km/s) and consistent with Galactic rotation (McGinn et al. 1989, Sellgren et al. 1990, Haller et al. 1996, Genzel et al. 1996). 

Eckart and Genzel (1996) have argued that the Hel stars also show a coherent pattern in their proper motions. Such a 
pattern is now confidently detected in the data (Fig. (|l^)). It is the origin for much of the tangential anisotropy discussed 
above. 11 of the 13 proper motion vectors for the emission line stars display a clockwise pattern, with only IRS 16 NE and 
IRS 16NW moving counter-clockwise. A number of authors have argued that most of the members of the IRS 16 complex 
(located between SgrA* and 4" east of it, and between 3.5" south of SgrA* and 1.5" north of it) belong to the early type 
cluster, with the Hel stars just being a sub-sample of the brightest emission line objects. This assertion is confirmed as well. 
Most of the brighter stars in the IRS 16 complex (m^ <13) show a clock- wise streaming pattern (Fig. (^) middle panel). 
In Fig. ( p^ (bottom panel) we overlay the proper motions vectors of Table 1 on the 0.05" resolution K-band Ghez et al. 
(1998) image of the SgrA* cluster. The preference of stars to be on radial/highly elliptical orbits that was discussed in the last 
section can be checked here from a graphical representation. Fig. (^sj) also suggests that the majority of stars in the SgrA* 
cluster have a similar projection of angular momentum along the line-of-sight, Lz, as the much brighter Hel stars and the 
IRS 16 cluster members. Speckle spectrophotometry (Genzel et al. 1997) and very recent high resolution VLT spectroscopy 
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Figure 13. Proper motion vectors on K-band images of the central star cluster. North is up and East is to the left. Top: proper motion 
vectors of all Hel stars with three velocities, overlayed on a grey scale SHARP map at 0.15" resolution. The length and direction of each 
arrow denotes the magnitude and direction of the proper motion for each star. For comparison a 500 km/s motion (for _Rq = S.Okpc) is 
shown in the upper left. To indicate the uncertainties in direction for the best stars, a shaded cone is shown for IRS 16C. The white cross 
denotes the location of SgrA*. The middle panel shows in addition the vectors for all stars in the IRS 16 complex that have K-magnitude 
< 13. This criterion likely selects primarily early type stars. Proper motion vectors within ^ 0.7" of SgrA* (the 'SgrA* cluster'), overlayed 
on the 0.05" resolution K-band map of Ghez et al. (1998) are given in the bottom panel. The asterisk marks the position of SgrA*. The 
la positional uncertainty of the radio source on the infrared map is between ±10 and ±30 milli-arseconds (Icr). The identifications of 
the sources (Sl-Sll, Table 1) are given. The length of the proper motion of SI corresponds to 1470 km/s for Rq = S.Okpc. 



© 0000 RAS, MNRAS 000, 000-000 



24 R. Genzel, C.Pichon, A.Eckart, O.Gerhard and T.Ott 



Table 3. Correction factor Mj^m/Mq for LM mass estimator versus /3 (horizontally) and s (vertically). Note the steep rise near /3 = 1 
for shallower density profiles. The function (AfLM/A^o)(/3, s) is illustrated in Fig. (tl4) 
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0.89 
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(Eckart et al. 1999) show that the brighter members of the SgrA* cluster lack 2.3-2.5/i m CO overtone absorption features 
and thus are clearly not late type stars. They are probably early type stars. If they are on the main sequence they would be 
of type BO to B2. The SgrA* cluster members thus are probably part of the early type star cluster but are on plunging, radial 
or very elliptical orbits. The only alternative explanation of the observed radial anisotropies and net angular momentum is 
that the SgrA* cluster stars are rotating as a group, like the Hel stars, but with a rotation axis lying in or near the plane of 
the sky. While this explanation seems relatively implausible, measurements of proper motion curvature and radial velocities 
are required to make a decisive test. 

However, the early type cluster cannot simply be modeled as an inclined, rotating thin disk. The fit of the best Keplerian 
disk model (inclination 40° , Vrot=200 r^""') to the Hel star velocities is poor. There are no Hel stars seen at Galactocentric 
radii greater than " 12". Because the Hel stars should be phase- mixed along their orbits (§2.2), a better description of their 
distribution, and perhaps the entire early type cluster, probably is a dynamically hot and geometrically thick, rotating torus 
at radii from 1" to 10" (0.039 to 0.39 pc). 

Most of the stars in the torus will have a fairly large angular momentum L and approximately the same sign of Lz. In the 
distribution of different L's there is a small fraction of stars, however, with much smaller L and still the same sign(Lz). This 
sub-population is necessarily small and may thus not contain very massive stars. The low-L stars are able to pass much closer 
to SgrA* than the majority of the early type star cluster. In our interpretation, it is these stars on largely radial, plunging 
orbits that make up the SgrA* cluster. As the bright, more massive stars are on average at larger distances from SgrA*, it is 
possible to detect fairly easily this central sub-sample of fainter, fast moving stars. One would expect to find the same types 
of stars also at larger radii from SgrA*. However, the present proper motion data sets are biased against such fainter stars 
because of the presence of the brighter early type stars (especially the IRS 16 cluster) and of late type stars at yet larger true 
radii. 

In summary of this section, we conclude that the majority of the Hel emission line stars and the bright (early type) 
stars in the IRS 16 cluster show a coherent clockwise and counter-Galactic rotation. Their circular (tangential) velocities 
dominate over their radial velocities. The young stars are arranged in a thick torus of mean radius «0.2 pc. This torus was 
presumably first formed ~7 to 9 million years ago when one or several infalling, tidally disrupted clouds collided and were 
highly compressed. This lead to an episode of active star formation in the central parsec. From the presence of bright AGB 
stars in the same region (Krabbe et al. 1995, Genzel et al. 1996, Blum, Sellgren and dePoy 1996) it is likely that there were 
other such phases of active star formation in the more distant past (a few lO'^ million years ago). 



6 PROJECTED MASS ESTIMATOR AND ANISOTROPY. 

Leonard and Merritt (1989) have shown that an anisotropy-independent, projected mass estimator can be constructed from 
radially complete proper motion data. Starting from the Jeans equation for a spherically symmetric, non-rotating system, 

GM{r) = -r-cr,^{dlogn(r)/dlogr + dlogCT,VdlogT + 2(1 - ol/al)} , (24) 

one can construct the spatially averaged, stellar tracer density (n{r)) weighted expectation value G{Ad{r)). This estimator, 
henceforth referred to as the Leonard-Merritt (LM) estimator, is independent of any assumptions about anisotropy, 

G{M{r))LM = {r{2a^r + 2cr|)) = (16/37r){2{i?a|) -f- {Ral)) . (25) 

Table 2 lists the LM-estimator obtained from 95 proper motion stars within 5" of SgrA*. The estimated mass is 2.9 ± 
0.4 X 10^ Mq. It is quite insensitive to the weighting scheme of the data. We also list LM-estimates for different projected 
annuli (although this is formally not appropriate, see below). 
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Figure 14. Leonard-Merritt mass estimates (in units of tiie true underlying mass) as a function of (5 for different power law slopes of 
clusters as labelled (the thicker curve corresponds topocr"^*). Note that for all power law clusters but that corresponding to p oc r~^, 
the mass estimate is typically biased and gives an overestimate or underestimate of the real mass. In particular note that the Bahcall 
Tremaine estimate (corresponding to /3 = 0) is also offset for all slopes but 3 and 1.62. For instance for /3 = and p (x r~^'^, the relative 
mass discrepency reads Mi^m / ^lo ~ 1-05 but for /3 = —5 it reaches 0.89 while a purely radial cluster would lead to an overestimation 
of 60 % as shown in Table 3. In practice our estimate of the mass of the black hole should be rescaled downward by about 5-10 % to 
account for this bias in the relevant radial range. 



For comparison, the Bahcall- Tremaine (1981, BT) estimator for an isotropic cluster around a point mass gives a mass of 
3.1 ±0.3 X 10® Mq for the 2x95 proper motions within R = 5" (Table 2). For purely radial orbits the mass would be twice as 
large. Note however that formally the BT estimator is defined for radial velocities only and as such the application to proper 
motions is inappropriate. The virial theorem mass estimate (VT) of the same data gives 2.5 — 2.6 x lO'' Mq (Table 2, see 
Bahcall and Tremaine 1981, or Genzel et al. 1996 for a discussion). The BT estimator requires prior knowledge of the orbit 
structure. In the region outside 3" from SgrA* where from our proper motion analysis the orbit structure is approximately 
isotropic, the agreement between all three estimators is fair (at somewhat more than la). 



6.1 Correction for the Leonard-Merritt mass estimate 

Unfortunately, the Leonard-Merritt mass estimate assumes that the cluster is of finite mass and that we have access to the 
full radial extent of the cluster. Here the density profile behaves roughly as a power law over the finite range of radii for which 
data is available. For such a mass model in a Keplerian potential the implementation of Leonard-Merritt mass estimate on 
concentric rings yields a biased (systematically offset) measure of the mass. Indeed the derivation of this estimator involves 
an integration by part of integrals of the form 



idpa (r) 
r , dr 



dr 



p(r)r-'V^(r)dr + [p(T')?"''a^(r)] 



(26) 



For finite mass systems the second term of the right hand side of Eq. 
still compute 



(l2q) vanishes. In the context of the Galactic centre we 
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I"r2 /■'■2 

/ p(r)r^(T^(r)dr/ / p{r)r^ a'^ {r)dr (27) 

J n ^ r\ 

over a finite radial range r\, T2 and even though both numerator and denominator diverge as ri ^ and ri oo, the ratio 
is well defined and equals the value corresponding to finite ri and r2- On the other hand 

[p{rya\r)]l\/ r p{rya\r)AT (28) 

is also finite but non zero except for s — —3. As a consequence the ratio Mlm/Mq will typically be a function of (3 and s, the 
slope of the local power law corresponding to the range for which data is available. A straightfoward calculation yields 

Mlm ^ 2^+- r{~s/2) ((3 ^s){l~s-P) Tj-s) - r(2 - s)) 

Mo 3 0F(l-s-2/3) r(-s/2-l/2)r(l/2-s/2)r(3/2-s/2) ^ ' 

Relevant relative mass estimates versus f3 for different power law index are shown in Fig. (^^. In practice Eq. (|2^) is used to 
correct for the offset in the measured Mlm- Table 3 gives a few values relevant for the Galactic centre. 
Note that Mlm = Mo for P — when s is —3 or the root of 

16r(s/2) ((1 -s) (3- s) r(s) -r(2 -s)) = 3 2'-' ^/^r(-s/2- l/2)r(3/2 - s/2f (30) 

which yields s « —1.62. More generally there is a non trivial curve (i.e. which differs from s — —3) corresponding to Mlm/A/q ~ 
1 in the (/3, s) plane. 

We conclude that the LM-estimator is not independent of /3 or s when applied to truncated data set even though each 
shell yields the same mass estimate for a Keplerian potential. We do need to estimate /3 and s independently to correct for 
the offset. Since (3 varies with radius for the Galactic centre the correction will affect the mass profile. 

The mass estimators are derived by averaging over the entire star cluster, while the observed stars in the Galactic centre 
are presumably part of a more extended stellar system. To test for possible systematic effects, we have therefore also carried 
out Monte Carlo simulations for the LM-estimator. Again, as in Section 4.2, we have used a power law distribution of tracer 
stars with n oc r~^'^ as for the kinematically measured stars. The left panel of Fig. (^5|) shows the median and quartiles of the 
distribution of Mlm values derived for many star cluster realizations, as a function of sample size and for /3* = 0.5,0,-0.5, 
and —3.5. The true mass of the central black hole that dominates the potential of these clusters is M = 1. The right panel of 
Fig. (^5|) shows, for A'' = 800 stars, the effect of estimating the central mass from five concentric annuli aranged linearly as a 
function of radius. Note that the mass profile is indeed flat (within the statistical uncertainties) as expected for a Keplerian 
potential and offset by the amount predicted by Eq. (psj). 

These simulations suggest that applying the LM-estimator to a central sub- volume of the actual star cluster around the 
black hole gives the correct hole mass if the distribution of orbits is strongly tangential, independent of power law slope (all 
radial shells should then be independently sufficient) . For isotropic and radially anisotropic orbit distributions and power law 
slopes near -2 the LM-estimator gives somewhat biased (too high) values for the central mass. The value of M = 2.9 x 10*' M© 
derived for the central mass from all stars inside 5" (an approximately overall isotropic sample) thus is likely systematically 
high by about 5 - 10%. 



6.2 Estimate of the Sun-Galactic centre distance Rq 

The expectation values of the first moments of the projected velocity dispersions are related to each other through their 
mutual dependence on the intrinsic radial and tangential velocity dispersions. One can write 

{i?<7,2) = (l/3){i?a|) + (2/3){i?4) . (31) 

The z- velocity is determined directly through the Doppler shifts of the stars. The R- and T-velocities depend on the assumed 
Sun-Galactic centre distance Rq. For a spatially and kinematically spherical system it is therefore possible to derive the 
distance to the Galactic centre from equation (31), without any prior assumptions on the anisotropy. The relationship is 

(J?0/8kpc) = i{Ral)s/{l/3{Ral)s + 2/3{Ra^)s}f'' . (32) 

Here ( )8 refers to the values calculated under the assumption that the Galactic centre distance is 8.0 kpc, as assumed 
for the proper motions in Table 1. Taking only those 32 stars for which we have all three velocity components we find 
Rq = 8.95 ± 1.6kpc. Taking all 104 proper motion stars within R < 8.8" and all 71 stars with z-velocities within the same 
projected radius we find Rq=8.2±0.9 kpc. The specific moment analysis in equation (32) as applied to these samples is 
appropriate if the motions are completely dominated by a central point mass. In that case Ra^ ~const and data points at 
different R (but the same quality) are appropriately given the same weight. In the Galactic centre the mass distribution is a 
sum of a central point mass and a near-isothermal stellar cluster of velocity dispersion (to=50 to 55 km/s derived from the 
stellar velocities outside the sphere of influence of the black hole (Genzel et al. 1996). It may thus be more appropriate to 
subtract ctq before computing the expectation values in equations (26) and (27). In that case we obtain Rq = 7.9 ± 0.85kpc. 
The difference between these two last estimates arises since the line-of-sight velocity data are biased to a larger (R) than the 
proper motions so that the effect of removing cro has a larger impact on the z-velocities. This differentially decreases slightly 
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Mass Profile 

AM/M versus N ,_, , , 




Figure 15. Left: Median and quartiles for the distribution of values for the mass estimator M^m (Eq. (|lO|)) from simulated proper 
motion samples, as a function of sample size and for several underlying anisotropy values = 0.5, 0, —0.5, —3.5 of the simulated star 
cluster. The true central mass is unity and the offset induced by the LM estimation is clearly visible there and in accordance with the 
theoretical prediction from Table {^. Mass determinations reliable at the 10% level require samples of at least 40, 35, 25, and 15 stars for 
/3* = 0.5, 0, —0.5, —3.5. Right: Central mass estimates Mlm derived in five annuli on the sky for a total number of 800 stars, anisotropy 
values are = 0.5, —0.5, and —3.5, as labelled (plain, dash-dot, long dash, dash) corresponding to (triangle, star, diamond, square). 



the distance estimate relative to that obtained for cro = 0. All errors do not contain a possible systematic term from deviations 
from spherical symmetry. 

Our analysis is in excellent agreement with other recent estimates for the Galactic centre distance which range between 
7.2 and 9.0 kpc with a best weighted average of 8.0±0.5 kpc (see the review of Reid 1993). The statistical uncertainty of our 
estimate rivals the best other methods available for determining Rq: cluster parallaxes through H2O maser proper motions, 
global modeling of the Galaxy, globular cluster dynamics, RR-Lyrae stars, Cepheids, planetary nebulae and OB stars in HII 
regions (Reid 1993) and clump giant stars (Paczynski & Stanek 1998). 



7 JEANS MODELING OF THE CENTRAL MASS DISTRIBUTION 

We have also carried out a full Jeans modeling of the data set, explicitly allowing for the anisotropy term in equation (p^. 
It is clear that the number of stars is still too small to unambiguously determine the radial profiles of anisotropy and mass 
for all different stellar components (and including rotation). Here we only give a simplified overall model which is consistent 
with all the data. Our model proceeds from a parameterized Ansatz for the different quantities, as described earlier in Genzel 
et al. (1996), 

n(r) = 2no/(7rn){l + (r-/ro)""}), (r) = ^^(r/ro)"^"'"' + (Tq and o-?(r-) = (t|o(7'Ao)"^"'" + o-q • (33) 

To compare to the observed surface density distribution E(i?), and observed velocity dispersions az{R), or{R) and (tt{R) 
the expressions in equation (P3|) were numerically integrated along the line of sight and weighted with the density distribution, 
as described in equations (jspand (^ . The data were averaged in annuli centered around SgrA* to yield 13 values for S(i?) 
between R=0.5 and 114", 8 values for (Jz{R) between R=2 and 100" and 5 values each for aR{R) and ariR) between R=0.4 
and 7.3". Best fit values for the 8 parameters in the expressions above were then determined from a minimization. They are 
listed in Table 5. The surface density measurements come from number counts with the SHARP speckle camera to m(K) = 15 
and are corrected for crowding and incompleteness (Schmitt 1995). The data points and their statistical errors are listed in 
Table 4. The R- and T- velocity dispersion values are from Table 2. The line of sight velocity dispersions are listed in Table 4 
as well. The two data points at R=2.2 and 8.5" are derived from the Hel star velocities in Table 1. In addition we have taken 
late type star velocity dispersions from Genzel et al. (1996) and references therein. The 13 surface density measurements in 
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Table 4 

Surface Densities and Line of Sight Velocity Dispersions 

projected distance from SgrA* m(K)<15 source surface density E 2AE/E 
(arcsecs) sources per square arcsec 



0.5 


3.5 


0.3 


1.5 


2.4 


0.3 


2.5 


2.2 


0.25 


4 


1.55 


0.3 


6 


1.06 


0.2 


8.5 


1 


0.2 


11.5 


0.57 


0.2 


16.5 


0.57 


0.3 


23 


0.44 


0.3 


28.8 


0.27 


0.32 


34.1 


0.22 


0.4 


73 


0.14 


0.4 


114.1 


0.08 


0.6 



projected distance from SgrA* 2 A(Oz)/Oz stellar type 

(arcsecs) (km/s) (km/s) 



2.2 


195 


0.35 


Hel 


8.5 


164 


0.9(a) 


Hel 


8.7 


102 


0.15 


late type 


10 


99 


0.21 


late type 


17.3 


72 


0.15 


late type 


20 


85 


0.35 


late type 


32 


68 


0.38 


late type 


100 


54 


0.22 


late type 



(a) given low weight to de-emphasize the dominant rotation signal of the early type stars at this radius 



Table 4 constrain the 3 parameters of the density distribution given in expression Eq. ( P3| ) very well. Likewise the 18 velocity 
dispersion measurements also give good constraints on the 5 parameters of the dispersion expressions in Eq. (^). 

Fig. shows the surface density and velocity data, along with the best fit model whose parameters are given in Table 5. 
Table 5 also lists the best anisotropic mass model, along with the logarithmic gradients and the 13 -anisotropy parameter that 
were used in the Jeans equation (Eq. (M )) to derive that mass model. 



8 DISCUSSION AND CONCLUSIONS 

The connected black crosses in Fig. ( |l7| ) depict the mass distribution obtained from the Jeans-model with anisotropy (and 
its (Icr) uncertainty). For comparison we also show the LM-mass estimators from Table 2, BT estimators for the z- and 
proper motions (this paper, Genzel et al. 1996, 1997, Ghez et al. 1998), the isotropic Jeans mass model of Genzel et al.(1996) 



© 0000 RAS, MNRAS 000, 000-000 



Stellar Dynamics in the Galactic centre: Pro])er Motions and Anisotropy 29 




offset from SgrA* (arcsecs) 

Figure 16. Surface density and velocity dispersions as a function of projected separation from SgrA*, along with the best fitting 
anisotropic Jeans model described in section 7. The observed m(K)<15 surface density counts and their Icr uncertainties are shown in 
the left inset (see also Table 4). These counts come from Schmitt (1995) and have been corrected for the effects of crowding and bright 
stars. The continuous curve is the model of Eq. (js^ with the parameters in Table 5, integrated along the line of sight as described in 
equations ( ^) and (^. The right inset shows observations of the projected tangential velocity dispersionss (T: triangles, this paper, 
Table 2), projected radial velocity dispersions (R: circles, this paper. Table 2) and line of sight velocity dispersions (z: rectangles with 
crosses. Table 4). The data for the line of sight dispersions come from this paper (Hel stars), from Genzel et al. (1996), Haller et al. 
(1996), Lindqvist et al. (1992), Sellgren et al. (1990) and McGinn et al. (1989). The best anisotropic Jeans model integrated along the 
line of sight (Table 5) gives the black, thin dashed curve for the projected tangential (T) data points, the red, thick continuous curve 
for the projected radial ( R ) data points and the blue, thin continuous curve for the line-of-sight (z) data points. The model assumes a 
distance of 8 kpc and gives a total of 21.6, or a reduced x^/N of 0.83. 
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Figure 17. Mass distribution in tlie central 10 pc of tlic Galaxy obtained from stellar and gas dynamics (for Ro = S.Okpc). The 
different data points (all with ilcr uncertainties) are as follows. Bold 'G' denotes mass estimates from the ionized and neutral gas 
dynamics (Serabyn and Lacy 1985, Glisten et al. 1987, Lacy et al. 1991, Herbst et al. 1993, Roberts and Goss 1993). Rectangles with 
crosses and down-pointing triangles denote the isotropic mass modeling of Genzel et al. (1996, 1997), including Jeans modeling of stellar 
radial velocities (early and late type stars, filled down pointing triangles) and Bahcall-Tremaine estimators of the NTT proper motions 
until 1996 (open down pointing triangles). Open rectangles (with crosses) are Bahcall-Tremaine estimators of the line of sight velocity 
data only. Open up-pointing triangles are the Bahcall-Tremaine estimators of the 1995-1996 proper motion data of Ghez et al. (1998). 
The new anisotropy independent mass estimates from the present work are given as filled black rectangles (Lconard-Merritt estimators 
of the proper motions in (Table 1, Table 2) and as large black crosses connected by a continuous curve (Jeans model, Table 5). For 
comparison several model curves are shown. The green, thick dashed curve represents the mass model for the (visible) stellar cluster 
{M/L{2fim) = 2,rcoro = 0.38pc, p(r = 0) = 3.5 X IO^Mqpc "3). The light blue, continuous curve is the sum of this stellar cluster, 
plus a point mass of 3.0 X 10^ Mq. The red, dashed curve is the sum of the visible stellar cluster, plus an a = 5 Plummer model 
(p(r) = p(0)(l + (r-/ro)^)~"/^) of a dark cluster of central density 3.7 X IO^^A/qpc^^ j.^ = 0.0058pc. 



and a few of the mass estimates determined from the gas motions (Guesten et al. 1987, Serabyn and Lacy 1985, Lacy et 
al. 1991, Herbst et al. 1993, Roberts et al. 1993). The mass inside of the innermost bin (r — O.Olpc) of our best Jeans 
model is 3.25 x W^Mq. It is consistent with the Leonard-Merritt mass estimator of the entire proper motion sample inside 
5" ,Mlm ~ 2.9 ± 0.4 X 10^ Mq (Table 2, see also Section 6). When corrected for the bias discussed above that mass becomes 
about 2.6 — 2.8 x 10® Mq. Systematic effects and the method of modelling dominate the accuracy to which the central mass 
can be determined. In Fig. ( |l7| ) we plot a central mass of 3.0 x 10® Mq. This value is a compromise between the bias-corrected 
LM-estimate and the Jeans estimate. Its overall (systematic plus statistical) uncertainty is ±0.5 x 10® Mq. It is reassuring that 
the results of our simple anisotropic modeling and of previous isotropic models are in good agreement within the respective 
errors. Our results confirm that the mass distribution is flat between 0.01 and 0.5pc. 
Nonetheless, there are still significant uncertainties in this analysis. 

(i) The parametric form of the model fitted to the data is not unique. 

(ii) We have so far not distinguished between early and late type stars. Yet it is fairly clear that early and late type stars 
have different spatial distributions and kinematics (see Plate 1 in Genzel et al. 1996). In the Jeans analysis, we require the 
density distribution and kinematics of an equilibrium tracer population. The proper motions stars are heavily influenced by 
spatial selection biases; thus it is not appropriate to use their inferred number density distribution in the Jeans equation. 
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Table 5 



Best Jeans Mass Model including Anisotropy ( for 8.0 kpc) 



n(U) 




Or 


X '^^ 






3.4 


1.25" 


1.8 


0.83 






oXr(O)) 




o,(r(0)) 


Ctm 


Oo 




280km/s 


1 


353 km/s 


0.5 


50 km/s 




r(pc) 


dlog n/ dlog r 


d log o / d log r 


2p(r) 


M(M®) 


2AM/M 


1.02E-02 


-1.03E-01 


-2.00E+00 


1 .33E+00 


3.25E+06 


7.50E-01 


1 .06E-02 


-1.09E-01 


-2.00E+00 


1.30E+00 


3.26E+06 


7.00E-01 


1.10E-02 


-1.17E-01 


-2.00E+00 


1.27E+00 


3.28E+06 


6.50E-01 


1.15E-02 


-1.25E-01 


-2.00E+00 


1 .24E+00 


3.29E+06 


6.50E-01 


1.19E-02 


-1.33E-01 


-1.99E+00 


1 .21 E+00 


3.29E+06 


8.00E-01 


1 .24E-02 


-1 .43E-01 


-1 .99E+00 


1.18E+00 


3.30E+06 


6.00E-01 


1.29E-02 


-1.53E-01 


-2.00E+00 


1.15E+00 


3.33E+06 


5.50E-01 


1 .35E-02 


-1.64E-01 


-2.00E+00 


1.11 E+00 


3.34E+06 


5.50E-01 


1 .40E-02 


-1 .74E-01 


-1 .99E+00 


1.08E+00 


3.35E+06 


5.50E-01 


1.46E-02 


-1.85E-01 


-1.99E+00 


1.04E+00 


3.36E+06 


5.00E-01 


1 .52E-02 


-1.98E-01 


-1.99E+00 


9.99E-01 


3.38E+06 


5.00E-01 


1 .58E-02 


-2.11E-01 


-1 .99E+00 


9.58E-01 


3.39E+06 


5.00E-01 


1.64E-02 


-2.24E-01 


-1.99E+00 


9.16E-01 


3.41 E+06 


4.70E-01 


1.71E-02 


-2.40E-01 


-1.99E+00 


8.72E-01 


3.42E+06 


4.50E-01 


1 .78E-02 


-2.55E-01 


-1 .99E+00 


8.26E-01 


3.44E+06 


4.30E-01 


1.85E-02 


-2.70E-01 


-1.99E+00 


7.78E-01 


3.45E+06 


4.20E-01 


1 .93E-02 


-2.88E-01 


-1.99E+00 


7.28E-01 


3.46E+06 


4.10E-01 


2.01 E-02 


-3.06E-01 


-1 .99E+00 


6.76E-01 


3.48E+06 


4.00E-01 


2.09E-02 


-3.24E-01 


-1.99E+00 


6.22E-01 


3.49E+06 


4.00E-01 


2.17E-02 


-3.44E-01 


-1.99E+00 


5.86E-01 


3.50E+06 


3.90E-01 


2.26E-02 


-3.64E-01 


-1 .99E+00 


5.08E-01 


3.51 E+06 


3.80E-01 


2.36E-02 


-3.86E-01 


-1.98E+00 


4.48E-01 


3.53E+06 


3.70E-01 


2.45E-02 


-4.08E-01 


-1 98E+00 


3 85E-01 


3 54E+06 


3 60E-01 


2.55E-02 


-4.31 E-01 


-1 .98E+00 


3.19E-01 


3.55E+06 


3.50E-01 


2.66E-02 


-4.55E-01 


-1.98E+00 


2.51 E-01 


3.56E+06 


3.40E-01 


2.76E-02 


-4.80E-01 


-1.98E+00 


1.80E-01 


3.57E+06 


3.30E-01 


2.88E-02 


-5.06E-01 


-1 .98E+00 


1.07E-01 


3.58E+06 


3.20E-01 


2.99E-02 


-5.33E-01 


-1.98E+00 


3.04E-02 


3.59E+06 


3.10E-01 


3.12E-02 


-5.60E-01 


-1.97E+00 


-4.89E-02 


3.60E+06 


3.00E-01 


3.24E-02 


-5.88E-01 


-1 .97E+00 


-1.31 E-01 


3.61 E+06 


3.00E-01 


3.38E-02 


-6.17E-01 


-1 .97E+00 


-2.17E-01 


3.61 E+06 


2.90E-01 


3.51 E-02 


-6.46E-01 


-1.97E+00 


-3.06E-01 


3.62E+06 


2.80E-01 


3.66E-02 


-6.76E-01 


-1.96E+00 


-3.98E-01 


3.63E+06 


2.70E-01 


3.81 E-02 


-7.08E-01 


-1.96E+00 


-4.94E-01 


3.63E+06 


2.70E-01 


3.96E-02 


-7.39E-01 


-1.96E+00 


-5.93E-01 


3.64E+06 


2.60E-01 


4.12E-02 


-7.69E-01 


-1.95E+00 


-6.96E-01 


3.64E+06 


2.60E-01 


4.29E-02 


-8.01 E-01 


-1.95E+00 


-8.03E-01 


3.64E+06 


2.50E-01 
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Table 5 (cont'd) 



r(pc) 


dloj 


; n/ dlog r 


d log d loi 


4.47E-02 


-8 


.33E-01 


-1.95E+00 


4.65E-02 


-8.66E-01 


-1 .94E+00 


4.84E-02 


-8 


..99E-01 


-1.94E+00 


5.04E-02 


-9.31 E-01 


-1.93E+00 


5.24E-02 


-9.63E-01 


-1 .93E+00 


5.46E-02 


-9.96E-01 


-1.92E+00 


5.68E-02 


-1 


.03E+00 


-1.91 E+00 


5.91 E-02 


-1 


.06E+00 


-1.91 E+00 


6.15E-02 


-1 


.09E+00 


-1.90E+00 


6.40E-02 


-1 


.12E+00 


-1.90E+00 


6.66E-02 


-1 


.15E+00 


-1 .89E+00 


6.93E-02 


-1 


.18E+00 


-1.88E+00 


7.22E-02 


-1 


.21E+00 


-1.87E+00 


7.51 E-02 


-1 


.24E+00 


-1 .86E+00 


7.82E-02 


-1 


.27E+00 


-1.85E+00 


8.14E-02 


-1 


.29E+00 


-1.84E+00 


8.47E-02 


-1 


.32E+00 


-1 .82E+00 


8.82E-02 


-1 


.34E+00 


-1.81 E+00 


9.18E-02 


-1 


.37E+00 


-1.79E+00 


9.55E-02 


-1 


.39E+00 


-1 .78E+00 


9.94E-02 


-1 


.41E+00 


-1.76E+00 


1.03E-01 


-1 


.43E+00 


-1.75E+00 


1 .08E-01 


-1 


.45E+00 


-1 .73E+00 


1.12E-01 


-1 


.47E+00 


-1.71 E+00 


1.17E-01 


-1 


.49E+00 


-1.69E+00 


1 .21 E-01 


-1 


.51 E+00 


-1 .67E+00 


1.26E-01 


-1 


.53E+00 


-1.64E+00 


1 .32E-01 


-1 


.54E+00 


-1.62E+00 


1 .37E-01 


-1 


.56E+00 


-1 .60E+00 


1.42E-01 


-1 


.57E+00 


-1.57E+00 


1.48E-01 


-1 


.59E+00 


-1.54E+00 


1 .54E-01 


-1 


.60E+00 


-1.51 E+00 


1.61 E-01 


-1 


.61 E+00 


-1.48E+00 


1 .67E-01 


-1 


.62E+00 


-1.45E+00 


1 .74E-01 


-1 


.64E+00 


-1 .42E+00 


1 .81 E-01 


-1 


.65E+00 


-1.38E+00 


1 .89E-01 


-1 


.66E+00 


-1.35E+00 


1 .96E-01 


-1 


.67E+00 


-1.31 E+00 


2.04E-01 


-1 


.67E+00 


-1 .28E+00 


2.13E-01 


-1 


.68E+00 


-1.24E+00 


2.21 E-01 


-1 


.69E+00 


-1.20E+00 


2.30E-01 


-1 


.70E+00 


-1.16E+00 


2.40E-01 


-1 


.70E+00 


-1.12E+00 


2.49E-01 


-1 


.71 E+00 


-1.08E+00 


2.60E-01 


-1 


.72E+00 


-1.04E+00 


2.70E-01 


-1 


.72E+00 


-1.00E+00 


2.81 E-01 


-1 


.73E+00 


-9.65E-01 


2.93E-01 


-1 


.73E+00 


-9.25E-01 



2|3(r) 


M(M®) 


2AM/M 


-9.13E-01 


3.64E+06 


2.50E-01 


-1.03E+00 


3.64E+06 


2.40E-01 


-1.15E+00 


3.64E+06 


2.40E-01 


-1.27E+00 


3.64E+06 


2.30E-01 


-1 .40E+00 


3.64E+06 


2.30E-01 


-1.53E+00 


3.63E+06 


2.20E-01 


-1.66E+00 


3.63E+06 


2.20E-01 


-1 .80E+00 


3.62E+06 


2.10E-01 


-1.95E+00 


3.62E+06 


2.10E-01 


-2.09E+00 


3.61 E+06 


2.10E-01 


-2.25E+00 


3.60E+06 


2.10E-01 


-2.40E+00 


3.60E+06 


2.00E-01 


-2.57E+00 


3.59E+06 


2.00E-01 


-2.73E+00 


3.58E+06 


2.00E-01 


-2.90E+00 


3.57E+06 


2.00E-01 


-3.08E+00 


3.56E+06 


2.00E-01 


-3.25E+00 


3.55E+06 


2.00E-01 


-3.44E+00 


3.54E+06 


2.00E-01 


-3.62E+00 


3.53E+06 


2.00E-01 


-3.81 E+00 


3.52E+06 


2.00E-01 


-4.00E+00 


3.51 E+06 


2.00E-01 


-4.20E+00 


3.50E+06 


2.00E-01 


-4.39E+00 


3.49E+06 


2.00E-01 


-4.59E+00 


3.48E+06 


2.00E-01 


-4.79E+00 


3.47E+06 


2.00E-01 


-4.99E+00 


3.46E+06 


2.00E-01 


-5.19E+00 


3.45E+06 


2.00E-01 


-5.39E+00 


3.44E+06 


2.00E-01 


-5.58E+00 


3.43E+06 


2.00E-01 


-5.78E+00 


3.42E+06 


2.00E-01 


-5.97E+00 


3.41 E+06 


2.00E-01 


-6.18E+00 


3.41 E+06 


2.00E-01 


-6.34E+00 


3.40E+06 


2.00E-01 


-6.52E+00 


3.39E+06 


2.00E-01 


-6.70E+00 


3.39E+06 


2.00E-01 


-6.86E+00 


3.38E+06 


2.00E-01 


-7.02E+00 


3.38E+06 


2.00E-01 


-7.16E+00 


3.37E+06 


2.10E-01 


-7.30E+00 


3.37E+06 


2.11 E-01 


-7.43E+00 


3.37E+06 


2.20E-01 


-7.54E+00 


3.36E+06 


2.30E-01 


-7.64E+00 


3.36E+06 


2.40E-01 


-7.74E+00 


3.36E+06 


2.50E-01 


-7.81 E+00 


3.36E+06 


2.60E-01 


-7.87E+00 


3.36E+06 


2.60E-01 


-7.92E+00 


3.36E+06 


2.70E-01 


-7.96E+00 


3.36E+06 


2.70E-01 


-7.98E+00 


3.37E+06 


2.80E-01 
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Table 5 (cont'd) 



r(pc) 


dlof 


; nj dlog r 


d log O^/ d log r 


2p(r) 


M(M®) 


2AM/M 


3.05E-01 


-1 


.74E+00 


-8.85E-01 


-7.98E+00 


3.37E+06 


2.80E-01 


3.17E-01 


-1 


.74E+00 


-8.46E-01 


-7.98E+00 


3.38E+06 


2.80E-01 


3.30E-01 


-1 


.74E+00 


-8.07E-01 


-7.95E+00 


3.38E+06 


2.90E-01 


3.43E-01 


-1 


.75E+00 


-7.69E-01 


-7.92E+00 


3.39E+06 


2.90E-01 


3.58E-01 


-1 


.75E+00 


-7.32E-01 


-7.87E+00 


3.39E+06 


3.10E-01 


3.72E-01 


-1 


.76E+00 


-6.95E-01 


-7.81 E+00 


3.40E+06 


3.30E-01 


3.87E-01 


-1 


.76E+00 


-6.59E-01 


-7.73E+00 


3.41 E+06 


3.50E-01 


4.03E-01 


-1 


.76E+00 


-6.24E-01 


-7.65E+00 


3.42E+06 


3.70E-01 


4.20E-01 


-1 


.76E+00 


-5.91 E-01 


-7.55E+00 


3.43E+06 


3.90E-01 


4.37E-01 


-1 


.77E+00 


-5.58E-01 


-7.45E+00 


3.44E+06 


4.10E-01 


4.54E-01 


-1 


.77E+00 


-5.26E-01 


-7.34E+00 


3.45E+06 


4.30E-01 


4.73E-01 


-1 


.77E+00 


-4.96E-01 


-7.21 E+00 


3.47E+06 


4.50E-01 


4.92E-01 


-1 


.77E+00 


-4.66E-01 


-7.08E+00 


3.48E+06 


5.00E-01 



Their role is to provide local velocity measurements for the population that they represent. The early-type stars contribute 
much to these kinematic measurements; if they are more centrally concentrated than the overall population measured by the 
SHARP number counts, this will have the effect of underestimating the central mass. 

(iii) The version of the Jeans equation we have used in expression neglects rotation. Yet we have discussed above the 
strong evidence for coherent motion of the Hel star cluster. The late type stars have only a small overall rotation. Including 
rotation and distinguishing in the analysis between late and early type stars would thus be desirable (as in Genzel et al. 
1996 for the isotropic case). Unfortunately this is not possible, because of the large number of free parameters (3 more for 
density distribution, approximately 8 more for velocity distribution) and the relatively poor constraints on a number of the 
parameters. There are no early type stars outside 11", there are very few late type stars inside 5" and the accuracy of the 
proper motion R- and T-velocity dispersions is low if all proper motions without a stellar type identification are discarded. We 
have run models with explicit inclusion of rotation but found it to be overall a poorer fit than the models without rotation. 
To deliberately deemphasize the rotation signature of the Hel star cluster we have arbitrarily given the z- velocity dispersion 
value at R = 8.5" (Table 4) a low weight. 

The central dark mass concentration is most likely a point mass. Any configuration other than a point mass must have 
a central density of p(0) > 3.7 x 10^^MQpc~"^and a core radius of ro < 5.8 milli-parsec . For this estimate we have adopted a 
Plummer model with a density profile that decreases as outside of the core radius. In a configuration with a point mass 
and the visible stellar cluster (p*(0) = 3.5 x 10^MQpc~^,a = 1.8, ro — 0.17pc) as the two main components of the mass 
distribution any additional mass within « 0.2 — 0.5pc of SgrA* must be less than « 1 x 10® Mq, or 32% of the point mass. If 
one takes the LM-mass distribution instead (Table 2), that limit would be between 1.1 and 2.2 x lO'' Mq. Backer (1996) has 
shown that the proper motion of SgrA* itself is < 16km s~^, or fifty to one hundred times smaller than the fast moving stars 
in its vicinity. Thus the mass enclosed within the radio size of SgrA* (r < lAU) is > 10^ or > lO'^ Mq, depending on whether 
the radio source is in momentum or energy equilibrium with the fast moving stars (Genzel et al. 1997, Reid 1999). Even the 
more conservative of these two limits implies a central density in excess of IO^^Mqpc-^ 

Our results confirm and strengthen recent work on the central mass distribution (cf. Eckart and Genzel 1996, 1997, Genzel 
et al. 1997, Genzel and Eckart 1997 , Ghez et al. 1998). From these papers and from Maoz (1998) it appears that the most 
likely configuration of the central mass concentration is a massive, but currently inactive black hole. With the parameters 
given above any dark cluster of stellar remnants (neutron stars, stellar black holes), low luminosity stars (e.g. white dwarfs) or 
sub-stellar objects would have a lifetime less than ~ 10^ years. This is much smaller than the ages of most of the stars in the 
Galactic centre, requiring that we happen to observe the Galactic centre in a highly improbable, special period. In addition, 
the very steep outer density distribution of such a dark cluster implied by the mass distribution in Fig. (0) (p « r-"with 
a > 5) is inconsistent with any known observed dynamical system. It is also inconsistent with the results of physical models, 
including those of core- collapsed clusters (see the discussion of Genzel et al. 1997). Maoz (1998) points out that the only 
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possible - albeit highly implausible - alternatives to a central black hole are a concentration of heavy bosons and a compact 
cluster of Ught (< O.OOSMq) 'mini'-black holes. 
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